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ABSTRACT 


This report contains the results of a fifteen month analytical and 
experimental study of the leakage rate of the pressurant gases (N^, he) 
and the propellant vapors (N 2 0^, N^) through bladder structures con- 
sisting of two layers of Teflon separated by a metallic foil diffusion 
barrier containing microscopic or larger holes* 

Results were obtained for the steady state leakage rate through 
circular holes and long rectangular openings in the barrier for arbitrary 
thicknesses of the two Teflon layers. The effect of hole shape and 
relative hole position on the leakage rate were studied. The transient 
problem was analyzed and it was shown that steady state calculations 
are adequate for estimating the leakage rate. A computer program 
entitled "Diffusion Analyzer Program" was developed to calculate the 
leakage rate, both transient and steady state. Finally, the analytical 
results were compared to experimentally determined values of the leakage 
rate through a model laminated bladder structure. The results of the 
analysis are in good agreement with experiment. 

The experimental effort (Part II of the Bladder Permeation Program) 
measured the solubility, diffusion coefficient and permeability of helium, 
nitrogen and nitrogen tetroxide vapor through Teflon TFE and FEP membranes. 
Data were obtained in the temperature range of 25 to 100°C at pressures 
ranging from near vacuum to about 20 atmospheres. Results of the experi- 
mental effort were compared with the limited data previously reported. 

As a verification to the applicability of results to actual bladder sys- 
tems, counter diffusion tests were performed with a laminated sample 
containing aluminum foil with a selected group of holes. 
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1.0 


INTRODUCTION 


1.1 Description of the Problem 

The objective of this program was to study the leakage of pressurant 
gases and propellant vapors through bladder structures composed of 
laminates of Teflon and thin sheets of stainless steel or aluminum foil 
serving as a diffusion barrier. As a result of pin-holes developed 
during manufacture and cracks which appear due to failure during flexure, 
the barriers are never perfect. If the barrier had no openings there 
would not be any leakage since the common pressurant gases as well as 
propellant vapors are not soluble to any degree in aluminum or stainless 
steel. Consequently, the purpose of the study can be stated as the 
calculation of the leakage rate through laminated bladder structures 
having holes in the diffusion barrier. 

A typical such laminated bladder structure consists of one layer 
of Teflon exposed to the pressurant gas, another exposed to the propel- 
lant, with a sheet of aluminum or stainless steel about 1 mil thick 
placed between. When there is an opening in the barrier, pressurant 
gas dissolves in the Teflon layer on the gas side and diffuses through 
the Teflon to the opening. It streams through the opening, or diffuses 
through, depending on whether or not the opening is filled with Teflon 
and dissolves in the other layer of Teflon. Subsequently, it diffuses 
through this layer of Teflon to the liquid side. Assuming no formation 
of gas between the liquid and the Teflon, the gas leaving the Teflon is 
dissolved in the liquid propellant and is transported by diffusion to 
the far reaches of the propellant. The driving force for this process 
is the gradient of the chemical potential, or the Gibbs free energy 
of the gas in solution. It can usually be stated as being merely due 
to the existence of a concentration gradient. Simultaneously, and by 
the same process, the propellant diffuses through the bladder into the 
pressurant gas space. 

Obviously, the problem, if taken in all its complexity, is quite 
difficult. A complete solution to the problem would require that one 
consider diffusion in the gas space, in the first Teflon layer, and 
through more than one irregularly shaped and spaced opening in the 
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barrier, followed by diffusion through the other Teflon layer, and then 
finally through the liquid propellant. In addition one must consider 
that propellant in molecular form simultaneously diffuses in the opposite 
direction, and that in general, the two diffusion processes are coupled, 
that is, the mass flux vector for the flow of the gas through the Teflon 
depends not only on the gradient of the concentration of the gas in 
Teflon but on that of the propellant as well. It seems reasonable to 
assume, however, that adequate engi^'ering results can be obtaine 1 using 
a somewhat simplified model. The model adopted is discussed in the 
following sections along with its justification. 


1,2 Multicomponent Counter Diffusion in Teflon 

The flow of one chemical component (the solute) In dilute solution 
in a second component (the solvent) obeys Ficks 1 Law^V as is well known. 
According to this law, the mass flux vector of the solute or the net 
amount of solute crossing unit area in unit time is proportional to the 
gradient of the concentration of the solute in the solution composed of 
the mixture of solute and solvent. Thus 


j = -pDVc 


( 1 - 2 . 1 ) 


where j denotes the mass flux vector, p the density of the solution, 

c the concentration of the solute in the solution, and D is the diffusion 

coefficient. The units of concentration chosen are those of mass 

3 

fraction, i.e., c = P./p, where p is the grams/cm cf species I present 
i i 3 1 

and p is the total gms/cm of the mixture, p « Ep^. Combining 
eq. (1.2.1) with the law of conservation of the mass of the solute 


"it * 


( 1 - 2 . 2 ) 


leads to the diffusion equation for a two component mixture 



V*pDVc 


( 1 . 2 . 3 ) 
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Equation (1.2.1) is valid for dilute solutions, in the absence of 
temperature gradients, pressure gradients, and chemical reactions. Thus 
eq. (1.2.3) is subject to the same restrictions. In fact both tempera- 
ture gradient diffusion and pressure gradient diffusion are very small 
effects so that for practical purposes eq. (1.2.3) is valid in the 
presence of quite large temperature and pressure gradients. The 
presence of a chemical reaction leads to a source term on the right side 
of eq. (1.2.3). 

When there are more than two components present, eq. (1.2.1) is 

more complex. The flow of component i depends on the gradient of the 

concentration of the other components present as well as that of com- 

. ( 2 ) 

ponent i 


• 4 P V C J 


(1.2.4) 


Thus in a mixture composed of pressurant gas with concentration c^, 
propellant with concentration c 2> and solid of concentration 1 - - Cj, 

we have 


— 

h 

-° D n vc i - 

pD 12 Vc 2 

(1.2.5) 


h = 

- pD 21 ,c l ‘ 

° D 22 Vc 2 

(1.2.6) 


where is the mass flux vector for pressurant gas, j 2 that for 
propellant. The law of conservation of mass then requires two equations 


= V-pDj^V^ + V*pD 12 Vc 2 


(1.2.7) 


p— = V- P D 21 V C;l + V-pD 22 Vc 2 


( 1 . 2 . 8 ) 


The coefficients D 12 and D 21 are called cross diffusion coefficients by 
some writers. The numbers D and D 22 are the direct diffusion coeffi- 
cients . 
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It is not possible to demonstrate theoretically the relative 
magnitudes of the cross diffusion and the direct diffusion effects. 

This can only be done by the actual measurement of all four coeffi- 
cients which would be an expensive and time consuming task. Instead, 
it will be assumed that the coupling between the two diffusive flows 
is negligible so that the concentration of each component obeys an 
equation of the form of (1.2.3). Thus we assume that the pressurant 
gas diffuses through the Teflon in one direction entirely independent 
of the presence of the pressurant diffusing in the opposite direction 
and conversely. Physically, one would expect a coupling to exist when 
the two components compete to a significant degree for the available 
lattice vacancies. Thus for dilute solutions, neglect of coupling 
effects should not lead to any serious error. 

This assumption can be tested experimentally by measuring the 

permeation rate for a given component in the absence of any counter- 

diffusing component and when one is present and comparing the two 

results. This was carried out as part of the experimental part of 

this program. It was found that the permeation rate of the pressurant 

gases was increased by the presence of counter diffusing » while that 

of N o 0. decreased. The direct diffusion coefficients can, however, be 
2 4 

adjusted to account for the effect to within the accuracy of the data. 

In the diffusion of through Teflon some chemical reaction 

apparently take place leading to swelling and discoloration of the 

Teflon. This effect could probably be accounted for by a direct 

diffusion coefficient which depends both on time and concentration of 

N„0, in Teflon. No effort was made under this contract to include 
2 4 

these effects in measured values of the diffusivity. However, the 
span of time over which D varies is probably small compared to times 
of interest for long term space missions so that sufficiently accurate 
calculations can be made using the asymptotic value of D(the diffusivity). 

Finally it should be noted that leakage through the Teflon could 
be due to leakage through capillaries in the Teflon, that is, Knudsen’s 
flow, in addition to diffusion. Such a transport process would^ not 
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obey the diffusion equation. Experiments carried out under this 
program, however, have found no evidence for the occurrence of such 
a process. (See Part II) 


1.3 Leakage Through Laminated Bladder Structures 

Suppose that a plane laminated bladder structure separates the 
gas and propellant spaces within a tank as shown in figure (1.1). 


Gas Space 


x = +h 


x = -b 



Liquid Space 

Figure (1.1) 

Let subscript 1 denote the diffusing component under consideration, 
for example, pressurant gas, leaking through the bladder into the 
propellant. As a result of the above assumptions concerning the 
diffusion of gas through the Teflon, we can write the following 
differential equations for the concentration c^ of the diffusing 
component . 

In the gas space: 


9c l 2 

at ' V c i 


h < x < 00 


(1.3.1) 


In the first Teflon layer: 


3c 2 

3T " D 1 V C 1 


0 < x < h 


(1.3.2) 


In the second Teflon layer: 


9c l 2 

3r - V C 1 


-b < x < 0 


(1.3.3) 
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In the liquid propellant: 


— = D V^c, - 00 < x < -b (1.3.4) 

8t 11 

Here subscript D refers to the diffusion coefficient of pressurant gas 
through a solution of pressurant gas and propellant molecules which 
have leaked through the bladder. is that for gas through the first 

layer of Teflon, D 2 for the second and D £ is that for pressurant gas 
through liquid propellant. The barrier, for convenience, is considered 
to have zero thickness. This assumption will not affect the accuracy 
of the calculations provided eqs. (1.3.2) and 1.3.3) are connected by 
the correct boundary conditions across the openings in the barrier. 

For finite gas and propellant spaces, the above equations have, 
strictly speaking, no non-trivial steady state solution. The pressurant 
gas, for example, leaks through the bladder structure into the liquid 
space, its concentration in the liquid space constantly increasing 
with time, until the entire liquid body becomes saturated with propel- 
lant gas. That is, the entire system reaches thermodynamic equilibrium, 
after which there is no longer any mass flux of either pressurant or 
propellant. However, for the nearly perfect barriers in which we are 
primarily interested, the openings are small and the leakage rate is also 
small. Quasi-steady state solutions are then possible which give a 
finite leak rate because the change in the concentration of the solute 
with time in either the pressurant or propellant space is very slow 
compared to the relaxation time of the bladder structure itself. 

Further, for small (microscopic) openings in the barrier, the region 
over which the distribution of concentration varies is very small 
compared with the volume of either the pressurant gas space or the 
liquid propellant space. These spaces may then, for purposes of 
calculation, be considered as semi-infinite space separated by a 
permeable bladder structure. 

In Part I - Appendix A, it is shown that under certain conditions 
the gas and liquid spaces are ignorable in solving eqs. (1.3.1 - 1.3.4). 
The total system is then replaced by one which considers only diffusion 
within the two Teflon layers, the concentration of solute being held at 
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a constant value on one surface and at zero value on the other. The 
conditions required are essentially that the permeation rates through 
the gas and liquid of the solute be much larger than that through the 
Teflon. Based on order of magnitude estimates of the diffusion parame- 
ters of gas through gas and gas through liquid and the values measured 
under this program for Teflon, the conditions seem to be well satisfied. 
The above system of equations are then replaced by 

0 < x < h (1. 3.5) 

b < x < 0 (1.3.6) 

where we now denote by c ^ the concentration of solute in Teflon layer 
two for convenience. The boundary conditions can now be chosen as 


9 c 


9c 2 J2. 

— - V C 2 


c^ = c Q - KjP 0 (x = h) (1.3.7) 

and 

c 2 = 0 (x = -b) (1.3.8) 

(3) 

and the conditions of thermodynamic equilibrium and conservation 
of diffusing particles across any opening in the barrier. is the 

solubility of solute in material 1 and Pq is the total system pressure. 
Outside the openings in the barrier the diffusive fluxes must vanish. 
The last condition is stated mathematically as 


9x 


9c, 


x=0 


9x 


- 0 


x=0 


(1.3.9) 


whenever (x,y) the coordinates in the plane lie outside the region of 
an opening In the barrier. Conservation of diffusing particles within 
the opening requires 


9c, 


P,D 


I"l8x 


9c, 


P 0 D, 


x=0 


2 29x 


x=0 


(1.3.10) 


whenever (x,y) lies inside the region of the opening. Here and 
refer to the total density of materials 1 and 2 respectively. 
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The condition of thermodynamic equilibrium at the interface between 

the two materials requires that the thermodynamic potential of the solute 

be equal at the interface. For dilute solutions of pressurant gas or 

(4) 

propellant in solids the thermodynamic potential takes the form 


M = RT2nc + f(p,T) (1.3.11) 

when V is the thermodynamic potential (Gibbs' Free Energy) , R is the 
gas constant, T the temperature and f(p,T) denotes a function of the 
pressure and temperature. From eq. (1.3.11) the condition of equilibrium 
is 


RT ^nc 1 + f 1 (p,T) = RT«-nc 2 + f 2 (p,T) 


or 



^(qcp.n-f^p.nj/RT . 


e- f l /RT 

e-V RI 


(1.3.12) 


The right side of eq, (1,3,12) is however just the ratio of the Henry 
Law solubility of component 1 in material 1 to that in material 2, 
Denoting the solubility by k we have 



as the condition of thermodynamic equilibrium at the interface between 
materials 1 and 2. 

Equations (1.3.7) through (1.3.13) indicate that the physical 
parameters needed to specify the leakage rate through a structure 
composed of two laminates of Teflon are the solubility and the 
diffusivity of the solute in each Teflon laminate. Since these 
quantities were not available in the literature, it was one of the 
objectives of this program to experimentally obtain these quantities as 
a function of pressure and temperature. The results of that phase of 
the program are reported in Part II. 
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2.0 STEADY STATE LEAKAGE THROUGH LAMINATED BLADDER STRUCTURES 

As was stated in section 1.3 steady state calculations of the 
leak rate through a laminated bladder structure are meaningful because 
for small leak rates the mean concentration of the diffusing species 
in either the pressurant or propellant tank changes significantly only 
over periods of the order of days whereas the transient period for a 
typical bladder structure is of the order of a few hours. Thus over 
long periods of time the concentration on the boundary of the bladder 
does not change appreciably so that the leakage rate can be calculated 
by assuming the concentration is constant. By assuming a maximum 
value of the concentration on one side, zero on the other, an upper 
bound for the steady state leak rate is obtained. Such results will be 
useful for engineering purposes and have the advantage of being con- 
servative, that is, they give the maximum leak rate. 


2.1 Leakage Rate Through Two Layers of Permeable Material 
With One Circular Hole in the Barrier 

This problem was worked out and reported in great detail in 
references 5 and 6 and will therefore only be outlined here. Consider 
Figure 2.1 which is a schematic of the cross-section of a laminated 
bladder structure containing a circular hole in the barrier foil 
of radius a. The concentration of the solute is held at c Q everywhere 
on the upper boundary, * 0 everywhere on the lower boundary. 

x 



c 

1 " °o 

x = +h 

1 

X'V'VX'V'V'V'V'V'V’V .... 

1 

Material 1 

r r rnnrpntrafinn r 



i 

VWWWWWW uat i J.&L 

Material 2 


I r=a Concentration c 

? X = -b 


Figure 2.1 
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As a consequence of the assumption of independent diffusion of the 
various species through the bladder material, the steady state concentra- 
tion of any given species obeys Laplace’s Equation. This follows from 
the usual diffusion equation by setting the time derivative equal to 
zero. If c^(r,x) is the concentration of the diffusing species in 
medium 1, c 2 (r,x) that in medium 2, then in cylindrical coordinates 

2 ft a 2 

c. -I 3 c.. 

— ir + — = 0 (0 < x < h) (2.1.1) 

Sr 2 r 9r 3x 2 - * 


, 3c I< 

$+±- 5 ^-+—^- 0 (-b < x < 0) (2.1.2) 

3 r 2 r 3r 3* 2 - - 

On the surface, x = h, the concentration c^(r,h) must go to a 
constant value, while on x * -b, c 2 (r,-b) is to vanish. Thus, 

c n (r,h) = c (all r) (2.1.3) 

1 * o 

c 2 (r,-b) = 0 (all r) (2.1.4) 

At the interface between the two materials there must be no diffusive 
flow except within the radius of the hole. Thus the presence of the 

=0 x = 0, r > a (2.1.5) 


=0 x = 0, r > a (2.1.6) 


barrier is simulated by 


3c. 


and 


3x 


3c„ 


3x 


and 


V 2 c, 
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Within the radius of the hole the conditions to be satisfied are 


3 

M 


w 


C;L (r,0) . ** Yc 2 (r,0) 


(x = 0, r < a) (2.1.7) 


9c (r,0) 

D l— 5— 


9c 1 (r,0) 

D 2— 5T~ 


(x = 0, r < a) (2.1.8)* 


Equations (2.1.7) and (2.1.8) state, respectively, that the concentra- 
tions in the adjoining materials within the radius of the hole are 
proportional, and that the number of particles is conserved (conservation 
of mass) . The proportionality of and c 2 requires the assumption of 
thermodynamic equilibrium at the interface between the two materials. 

The proportionality constant, Y, is known as Nemst's constant and is 
equal to the ratio <^/ < 2 of the solubilities of the solute in the two 
materials. It is not strictly constant but depends on temperature and 
pressure, principally the temperature^ \ Finally, note that c q in 
eq. (2.1.3) is not the concentration of diffusing component (i.e.. 

He gas) in the space above x = h, but is the equilibrium concentration 
of this substance in material 1 when it is exposed to a bath of pure 
diffusing component. Thus, c^ must be calculated from thermodynamic 
equilibrium data which will be obtained experimentally. 

For convenience, we solve not for c^(r,x) in region 1, but for the 
function <Kr,x) defined by 


Kr,x) = c q - c 1 (r,x) 


(2.1.9) 


This function has the advantage of vanishing at x = h, and enables one 
to avoid the use of non— convergent integrals. Obviously, it satisfies 
Laplace's equation. 


ar 2 r3r a * 2 


0 < x < h 


( 2 . 1 . 10 ) 


*We assume the two Teflon materials have the same density, = p 
The results are easily modified to include the more general case. 
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Taking the zero-order Handel Transform of Equation (2.1.10) removes the 
r dependence. The transform is defined as 


* (k,x) 


f°° 

- n 


’Kr.x) 


r J (kr) dr 
o 


( 2 . 1 . 11 ) 


Using eq. (2.1.11) in eq. (2.1.10) gives as the differential equation 
obeyed by ^(k,x) 

^ - k 2 ^(k,x) = 0 (2.1.12) 

dx 

The solution to this equation which vanishes at the boundary x = h is 

+< k **> - V k > IlnhkiP 1 (2 - 1 - 13) 


where A^(k) is an unknown function of k. Of course, the denominator of 
eq, (2,1.13) can be absorbed in the function A^(k) if desired. By the 
inversion theorem for Hankel transforms. 


ij>(r,x) = k A 1 (k) g - 1 ^h^kh ^' J 0 (kr)dk (2.1.14) 

rv 


so that 


c 1 (r,x) = c - 


f k A l (k) J o (kr)dk (2a - 15) 


In the same manner it is found that the solution for x < a is given by 


c 2 (r,x) 


r 

kA 2 (k) 

• n 


sinh k(b+x) 
sinh kb 


J (kr) dk 
o 


(2.1.16) 
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The two sets of boundary conditions (2.1.5), (2.1.6) and (2.1.8) 
can be combined to state that 


for all r. 
then gives 


3c 1 

W 


9c, 


D, 


2 9x 


x = 0 


Application of this condition to eqs. (2.1.15) 


(2.1.17) 
and (2.1.16) 


T> 1 [ k 2 A x (k) ctnh kh J q (kr)dk = D 2 | ctnh kb J q (kr)dk 

' o ■'o 

(2.1.18) 

Since this equation is valid for all r, the Hankel inversion theorem 
states that 


A 1 (k) ctnh kh = D 2 A 2 (k) ctnh kb 


or 


* ft,-) m ctnh kh /^\ 

V' D 2 ctnh kb V ' 


(2.1.19) 


for all k. Using now the condition, eq. (2.1.7) and eq. (2.1.19) 
we find that for r < a. 


r kA.(k) ( y D A 

— tanh kh + — — tanh kb J (kr)dk = 

tanh kh ^ D 2 Jo 


= c (2.1.20) 
o 


while for r > a, from eq. (2.1.5), 


k 2 A x (k) 


tanh kh o 


J (kr) dk = 0 


( 2 . 1 . 21 ) 
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The latter equation being required in order that the current vanish for 
r > a. Putting 


lc A^k) 

tanh kh 


= f(k) 


reduces eqs. ( 2 . 1 . 20 ) and ( 2 . 1 . 21 ) to the pair 


f(k) 

k 


Y 1 1 

tanh kh + tanh kb J Q (kr) dk = c q , (r < a) 

2 ( 2 . 1 . 22 ) 


f(k) J (kr) dk = 0 , 
o 


r > a 


(2.1.23) 


Substituting P = 7, k = 7, a = 7 and S = 7 we have the set 


f°° f( ) f YD. ) 

r ^ u - tanh otu + — — tanh f5u J (pu)du = 
_ u [ D 2 Jo 


c p < 1 
o 


(2.1.24) 


and 


f (u) J (pu)du = 0 
o 


p > 1 (2.1.25) 


These two equations are of the form 

G(u) f (u) J (pu)du = c 


p < 1 (2.1.26) 


f(u) J (pu)du = 0 
o 


p > 1 (2.1.27) 


where 


G(u) - t 


tanh cm + - — tanh 6 u 
2 


(2.1.28) 
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A general method of solution, due to Tranter , is available for this 
type of dual integral equations. It is based on the fact that if f(u) 

Is expanded as a series of Bessel functions of ascending order, 
eq. (2.1.27) is identically satisfied. The coefficients in the expansion 
have then to be chosen so that eq, (2.1.26) is satisfied. The procedure 
is complex, however, and can only be carried out to yield closed form 
solutions in a few cases. Putting 


f(u) - u 1 -" l a m W ») (2.1.29) 

m-o 


eq. (2.1.27) is satisfied due to the discontinuous nature of the 
resulting integrals. The substitution of eq. (2.1.29) into eq. (2.1.26) 
leads to an infinite set of linear inhomogeneous equations for the 
expansion coefficients a^. Tranter presents an iterative technique for 
solving these equations which works well in some cases. A slight 
modification of his technique works well in this case also so long as 
a and 6 are large. Large a and 8 corresponds to the case where the 
thickness of the permeable laminate is much greater than the hole 
radius. This is a practically important case. 

The details of the calculation of the coefficients a m and 
subsequently the leak rate are given in reference 6. Here we 
quote and discuss the result for the total flow rate W which is given 
by 
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(2.1.30) 


accurate to terms of third degree. 

The quantity W/4D^ac Q p depends on three parameters. These are 

the ratios a/h, a/b, and the factor yD^/D^. Figure 2.2 shows a typical 

set of curves for various values of y * b/h, the ratio of the thickness 

of the two slabs and for yD^/'D^ = 0.5. The parameter yD^/^ is, in 

fact, the ratio of the permeability of material 1 to material 2. A 

(9) 

more complete set of curves is available in the Design Guide 

The result is expected to be valid for ratios of hole diameter 
to laminate thickness as large as 0.5, perhaps even larger. The 
reason for this is that eq. (2.1.30) is not very sensitive to the terms 
nonlinear in a/h or a/b since the coefficient of these terms is small 
compared to that of the linear terms. The quantity c(3) in eq. (2.1.30) 
denotes the Riemann Zeta function of argument 3.0 and has the value 
1.202 nearly. 
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Figure 2.2 Dependence of Flow Rate on the Ratio of Hole Radius to 

Slab Thickness 

Nomenclature: 


hole radius 

Teflon thickness, side 1 
Teflon thickness, side 2 
ratio of Teflon thicknesses 
diffusion coefficients 
ratio solubility constants 


mass flow rate 
Teflon density 

Concentration of diffusion gas 
absorbed by Teflon at Teflon- 
gas space interface, dimensions 
of mass of gas per unit mass of 
Teflon 




To use the results to calculate W, the total flow rate per unit 

time, one must have the value of y, the two diffusion coefficients, the 

hole radius, and the thickness of the two slabs. In addition, one needs 

to know the value of the concentration c^, on one surface. If, for 

example, it is desired to calculate the leakage of Nitrogen pressurant 

gas through the bladder structure to the propellant side, the value of 

c is obtained from the Henry law constant and the pressure of the gas 
o 

by c q = K^p , where is the Henry law constant for material 1 and p 
is the pressure of the gas pressurant. The constant y is the ratio of 
the two Henry law constants for the two permeable materials, that is, 

y = k 1 / |< 2* 

Given y and the two diffusion coefficients one can estimate the 
parameter yD^/D 2 - From the known value of a/h, the ratio of hole radius 
to gas side slab thickness, and y = b/h, the ratio of the liquid side 
slab thickness to that of the gas side one, the quantity W/4D^ac Q p can 
be obtained from the curves. 

If the value of does not exactly match any of the plotted 

values, linear interpolation between two charts can be used. From the 
numerical value of W/4D^ac Q p thus obtained, one obtains the value of W. 
The value of p in W/4D^ac Q p refers to the density of the solvent, 
that is, the permeable material. It has been assumed that the density 
of the two layers is the same (eq. 2.1.8). If not, we replace by 
p^D^ and D 2 by P 2^2 ex P ress * on 1 + yD^/^ appears in 

eq. (2.1.30) and replace p by P 2 in W/4D^ac o p. 
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2.2 Leakage Rate Through Two Layers of Permeable Material 

With an Infinite Slit in the Barrier _ 

In the preceding section the problem of the leakage rate through 
a bladder structure composed of two laminates of permeable material 
separated by an impermeable barrier with one circular hole was solved. 
The solution obtained was rigorous provided the radius of the hole Is 
less than the thickness of either laminate. The only other simple 
geometry Is obtained when the opening is an infinite slit. Such an 
opening corresponds to the real case of a rectangular opening which 
Is very narrow so that the effects of the ends can be neglected. 

It has not been possible to find a rigorous solution to this 
problem in the manner used for circular openings in section 2.1, 
although several attempts were made. An approximate solution can be 
constructed however by assuming a form for the leakage current across 
the slit opening and equating the average concentration values on each 
side of the slit proportioned by the ratio of the solubilities. A 
good assumption for the form of the leakage current is that for 
leakage through a slit into a semi-infinite medium. By comparing the 
results of such a calculation with those obtained using the digital 
computer, we can assess the validity of the assumption. The analytical 
calculations have advantages over those obtained using the digital 
program since the results can be tabulated In terms of dimensionless 
parameters as for a circular hole and will be available without the 
expense and effort of a computer run for any given case. Further, the 
assumptions stated above are quite reasonable so that the results are 
not expected to differ significantly from those of a more rigorous 
solution. 


A cross section of the structure is shown in figure (2.3). The 

* y 


c.. = c 
1 o 


y = +h 


Barrier 




Material 1, Concentration c 1 


» » »» » » ? - » — ;<«UUUU 


"I |' 

x= -a x~ 


Material 2, Concentration c„ 
a , 2 

y - -b 


0 


Figure 2.3 
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width of the slit is 2a. In material 1, for y > 0, the concentration 
of the solute c^ satisfies Laplace’s equation in two dimensions 




( 2 . 2 . 1 ) 


and takes on the boundary condition c^ = at y = +h for all x. 
Specification of the flow through the slit leads to a complete solution 
for y >0. Along y = 0 we put 


3c 1 (x,0) 

p i D i — 37“ = w(x) ’ 


X < a 


= 0 


x > a 


( 2 . 2 . 2 ) 


As stated above, the form of w(x) for flow into or out of a semi - 

pi 2~ 

infinite medium through a slit of width 2a is w(x) « 1/ /a - x 
We determine the proportionality constant so that the integral of w(x) 
over the width of the slit is equal to the total steady state flow 
rate W (gms/sec per unit depth of slit)* Since 


r 

w 

J- a 


(x) dx 


f+a 

i-a 



ir 


this is accomplished by putting 

w(x) = - — 1 - (2.2.3) 

” -yfl 2 
Va - x 

The second boundary condition then reads 

3Ci<*, 0) „ j 

p i D i ay 7 n I l x l <a 

/ a - x 

(2.2.4) 

* 0 |x| > a 


2 - 12 


As in section 2.1, replace c^ by c^ - »Kx,y) . The conditions 
by \f)(x,y ) are 

2 2 

i_J? + ±J> = o 
2 2 
3x ay 


i(j(x,h) = 0 


aii>(x.o) _ w i 
• p l D l ay - , 


|x I < a 
lx I > a 


Since the problem is completely symmetrical about x = 0, 
Fourier cosine transform can be employed. This is defined as 


;^(k,y) = c 1 ( x »y) cos kx dx 

* r\ 


with inverse 


c i (x 'y ) ■ •; 


c^(k,y) cos kx dk 


The cosine transform of iK x ,y) satisfies 


d - k 2 #k,y) - 0 

dy 


with solution vanishing on y = h. 


,, \ . /i\ sinh k(h - y) 

>Kk,y) =* A^k) sinh fch 


© 


satisfied 

(2.2.5) 

(2.2.6) 

(2.2.7) 

the 


( 2 . 2 . 8 ) 

(2.2.9) 
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Taking the cosine transform of eq. (2.2,7) gives 


\\ \ cos kx dx = -? } 
o y c 

d^(k.O) 


cos kx dx 


The left hand side is P,D n - while the value of the right side 

w 1 1 dy , q \ 

is - — J (ka) where J (ka) is the zero order Bessel Function v . Using 
2 o o 

eq. (2.2.9) we obtain 


kA^(k) ctnh kh = ^ J o (ka) 


( 2 . 2 . 10 ) 


from which 


4»( k >y) 


W tanh kh * sinh k(h-y) 

k o^ a) sinh kh 


( 2 . 2 . 11 ) 


Taking the inverse transform gives c^(x,y) 

W 


c^x.y) - c o - Pi „ 


, , , sinh k(h-y) T x , dk 

tanh kh . - J ( ka ) cos kx T" 

sinh kh o k 

o 

( 2 . 2 . 12 ) 


The value of W cannot be determined from this result. It must be 
calculated by obtaining a similar result for y < 0 and then requiring 
that the concentrations be propotional within the slit opening. 

The solution for c 2 (k,y) which vanishes at y - - b is 


c 2 (k,y) 


= A 2 (k) 


sinh k(b+y) 
sinh kb 


(2.2.13) 


Again the condition given by eq. (2.2.4) must be satisfied which gives 


, ,, x W tanh kb , „ x 

V k) ' — — J o (ka) 


(2.2.14) 


and upon using the value of A^(k) and applying the inversion theorem 
we obtain 


W 


c 2^ x »y) - rrr 


2 2 


tanh kb s * n |? J (ka) cos kx ^ 

sinh kb o k 


(2.2.15) 
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The flow rate W is now determined by requiring that 


c^x.O) = yc 2 (x,0) 0 < x < a (2.2.16) 

as in section 2.1, where y is the ratio of the permeabilities of the 
two materials. From eqs. (2.2.12) and (2.2.15) we obtain 


c - 


W 


r°° 


* p i d i 


dk 


tanh kh J (ka) cos kx . 

o k 


(2.2.17) 


= T 


W 


r 


" p 2 D 2 


tanh kh J (ka) cos kx 4^ 0 < x < a 

o k - - 


Clearly this will lead to a unique determination of W only if 
the integrals in eqs. (2.2.17) are independent of x in the interval 
0 < x < a. That this is not the case can be determined from the value 
of the integrals for small h and b. In that case tanh kh can be 
replaced in the integral by kh and the value of the resulting integral 

n — 7 

is h/Va - x for 0 < x < a. The only way out of this difficulty is 
to average the equation over the slit opening by integrating out x. 
This gives 


W 


r 


c — 

o Trap..D 


11 


dk 


tanh kh J (ka) sin ka —7 

O 1 ^ 

o k 


= r 


[ tanh kh J (ka) 

wap 2 D 2 Jo 0 k 2 


(2.2.18) 


which is independent of x. Solving for W gives 


W = 


Tiap 1 D 1 c 
1 1 o 


(tanh kh + y tanh kh) J (ka) s ^ n ^k 
o 2 2 k 


(2.2.19) 
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The integrals are evaluated in Part 1, Appendix B by the Theory of 
Residues. The integrals are of the same type and depend only on either 
h and a or b and a. Using the result obtained in Appendix B, we find 

2ap 1 D 1 c 

w LJ-S-S 

o YP,D, 9 

h ( x - G ( h> ) + b C 1 -T G( t^ (2 - 2 - 20) 

tt n p 2 2 n 

as the leakage rate in mass per unit time per unit length through a 
bladder composed of two laminates of thickness h and b respectively 
with a slit of width 2a in the barrier. The function G is given by 


G(a) 



n=o 


I [ (n-Hr) ttu] 
o l 

(n+i ) 2 


- (n-hr) ira 
e l 


( 2 . 2 . 21 ) 


when I denotes the zero order modified Bessel function, 
o 

The quantity (Zap^D^c^) / [h-Kyp^D^^Dpb] is the leakage rate 
through a two laminate bladder of width 2a and unit depth in the absence 
of a barrier. In order for the expression (2.3.20) for the leakage rate 
to give the correct value for a >> h; a >> b it is necessary that G 
vanish for large values of the argument a/h; a/b. This is indeed the 
case (See Appendix B) . 

Figure 2.4 shows a typical curve of W/p^D^c q plotted as a function 
of the ratios a/h and y = b/h. A more complete collection of such curves 
will be found in reference 9. Note that for a given value of 
yp^D^/p^D^ and y the curves become linear in a, the slit half width, 
when a/h is of the order of 0.5. This indicates that when the total 
width of the slit, 2a, is of the order of the laminate thickness the 
effect of slit edge effects are negligible and the leakage rate can be 
calculated using the formula quoted above for leakage through a. strip 
composed of two permeable slabs without a barrier, the width of the 
strip, of course, being 2a. Since the results were obtained by an 

approximate method we have compared them with numerical calculations 
using the Diffusion Analyzer Program, DAP. Figure 2.5 gives a compari- 
son of the two results. It can be seen that the approximate analysis 
Is quite adequate. 
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2a/h 

Figure 2.4 Diffusion Rates Through Laminated Bladder with 
Infinite Slit in Barrier 
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Computer Results 

Analytical Results 

Note: See Figure 2.4 for 

definition of symbols 


2 1.6 2.0 2.4 


ter Results with Analytical 
r Through Laminated Bladder 
in Barrier 



2.3 Effective Area of Influence of Holes 


Calculations of the leakage rate through a single hole in the 
barrier can be used for multiple holes provided these holes are 
sufficiently far apart. The calculation of section 2.1 can be used 
to determine the distance by which circular holes must be separated 
in order that leakage through two holes can be considered independent. 

Referring to section (2.1), the concentration in the upper 
material (material 1) is given by 


c 1 (r,x) - c q - f(r,x) 
and from eq. (2.1.14) <J> is given by 


(2.3.1) 


^ (r ,x) 


r oo 

kA x (k) 

o 


sinh k(h-x) 
sinh kh 


J (kr) dk 
o 


(2.3.2) 


when the value of A^(k) is known from reference 6. If there were no 
hole in the barrier the concentration throughout material 1 would 
equal c . Therefore the radius at which >Kr,o) effectively vanishes 
is a good measure of the area effected by the hole. Inserting the 
value of A^(k) from reference 6 we have 


Kr,o) 


kA 1 (k) 


J (kh) 
o 


dk - 


a a 
o 




1/2 J tanh kh 

, 1/2 


J l/ 2 ( k a) 


J (kr)dk 
o 


+ 


a l a 


1/2 


tanh kh 
k 1/2 


J 5 ^ 2 (ka) J Q (kr) dk 


(2.3.3) 


Since a^ is much smaller than a^ and also the second integral is 
smaller than the first, we can write 


H r,o) 


a a 
o 


1/2 


■ CD 


J 

o 


tanh kh 
k 1/2 


J /0 (ka) J (kr) dk 
1/2 o 


(2.3.4) 
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closely. Since J, ,„(ka) = — 

1 / 2 y TT 

*(r,o) a o f - 

J r\ 


sin ka 


, eq. (2.3.4) is equivalent to 


— sin ka J (kr)dk 
k o 


Putting ka = u reduces eq. (2.3.5) to 


o* r 


tanh £u sin u J (yu) du 
o 


where 3 = h/a and y = r/a. Finally, using the expansion 


tanh 3u = 1 + 2 } ' (-1) S e 


gives 


: - o) -£% | ( i + 2 


(-1) S e ^ s ® u j gin u J Q (yu) (2.3.7) 


Integrating term by term after interchanging the order of integration 


and summation gives 


:r - o) -M%\ 

* c 


sin u J (yu) — 
o u 


+ J— a T] (-1) [ e sin u J (yu) — 

V TT O J O U 

s = l J O 


(2.3.8) 


We are interested in y > 1. In this case the value of the first term 
is (8) 

[2 . -l.a. (2 fa 1 a^ 3 a "* . 1 

\ — a sin (— ) = ./ — a ”+T“T + 777“T + • • * 

V tt o r V tt o J^r 6 ^3 40 ^5 J 
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The second integral is also known. 


-2s n u sin u T ^ 

e 3 f— J (yu) du 

u o 


sin 


•i f 

lVasV + (i + -) 2 +v4s 2 3 2 + (i - ~) 2 - 


We thus find 
*(r 


» o) = Vf a 0 [f + £ % + & % + • • •] 

u r r J 


* a ° s 5 


(-1) sin 


4 _ 1 1 

L^sV + (1 + -> 2 -p/^V + (i - f) 2 -* 


(2.3.9) 


The leading term in is independent of hole radius and ratio 
of slab thicknesses to hole radius. Therefore we calculate the ratio 
ip( r,o)/a Q as a function of a/r for various values of 3. Figure 2.5 
shows a plot of i|j(r,o)/a^ as a function of a/r from eq. (2.3.9). It 
can be seen that the value of this quantity at a/r = 5.0 is about 10% 
of its value at a/r = 1.0, the hole edge. This ratio is not strongly 
dependent on the value 3. Therefore, we conclude that 10-20 hole 
radii should be about the maximum extent over which the effect of a 
hole is felt. The same conclusion is, of course, valid for the 
second material. 
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r/a - Dimensionless Radius 








2 . 4 Leakage Through Periodic Slits 

There is one other case which can be treated analytically in an 

approximate manner and which is of interest for determining the effect 

of adjacent holes on the leakage through each. This is the case of a 

laminated bladder structure constructed of laminates of permeable 

material and a barrier with periodic slits. A cross section of such 

a structure is shown in figure (2.7). The concentration on the barrier 

side is held at c while that on the opposite side is zero. Each of the 
o 

slits is of equal width, 2a, and the intervening barrier is of width 
2(d -a). 

| 

I c = 0 

y = h 

Permeable Material 

"l-U ' 

Figure 2.7 

The concentration obeys Laplace T s Equation in two dimensions (x,y) . 



+ = 0 (2.4.1) 

3x 2 3y 

Since the concentration must be a periodic function with period 2d 
we expand c(x,y) as follows: 

oo 

c(x,y) = f n (y) (2.4.2) 

n=o 

where f (y) is an arbitrary function of y. Substitution of eq. (2.4.2) 
n nTTX 

in (2.4.1) and use of the linear independence of the set cos — gives 

as the equation of the f 

d 2 f - 

— ^ - (~r) 2 f (y) = ° (2.4.3) 

dy 2 d n 


■ Ti 
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k-nui- 


with solution vanishing on y = h of 


sinh ^j-(h-y) 

f n (y) A n mrh 

sinh — -j— 


(2.4.4) 


where the A are arbitrary constants. Hence 
n 


c(x,y) = 2 A n 


sinh ^r(h-y) 

d mrx 

-r- cos — j~ 

, , nirh d 


(2.4.5) 


The current flow is 


9c(x,o) 


PD 2^ T A n 


d ntx 

mrh COS d 


(2.4.6) 


The flow rate in the interval 0 < x < d can be expanded in a Fourier 
cosine series. If W is the total leakage rate through the slit 
-a < x < +a and if we assume the flow is uniform then W/2a is the local 
current (per unit depth of bladder). The current vanishes, of course, 
for a < x < d. 


00 

S ' E B n cos 


(2.4.7) 


where 


„ 1W.1W_ _W_ 

B o “ a j 2a dx d 2a 2d 
J o 


2 W nTrx 

d , 2a C ° S d 
4 o 


ittx , W 1 nira 

— : — dx = sin —j— 

d rra n d 


Therefore 


w 

W , W 1 . .nira. 

2d + « If n sln > cos 


(2.4.8) 
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Equating the two expressions for the current gives 


A 

c 

h ' d 


pD <-£ + 4 


E , , mTh % ,u iix x 

n A ctnh(— — ) cos (— -j-) 
n d d 


/ n tdc x 


n=l 


Thus 


W_ h 

o " 2pD d 


W . W V 1 . ^mra-v ,nTTX x 

T7 + IT" 2- — sin (—7-) cos (-j-) 

2d n a , n d d 

n=l 

. , ,mrh x 

w d sin (— ) tanh > 


n tt 2 pD 3 


n 


(2.4.9) 


giving at y = 0 


, ^ w n w a 

c(x>o) 2- 

M u p D n=l 


, /mra x ^ /niih^ 
00 sin (—5—) tanh (— 'j— ) 


n 


— cos ^ (2.4.9) 

a 


We now equate the average value of c(x,o) over 0 < x < a to c , the 
concentration there. This gives 

,2 

W h . W c 

c = 


W h , _W d_ 

“o ” 2pD d _3 ^ 2 


7T pD a 


, nira , mrh 
sin — tanh — n „ a 

3 Sl " T 


(2.4.11) 


Solving for W we obtain 


2pD c 


W = 


2 OO 

h + Ad _ y 

d 3 2 Z. 

L tt a n= 


, , ntra>.2 , 
(sin — t— ) tanh 

£ — Hr— 

n=l n 


,nrrh. " 
2 d ’ 


(2.4.12) 


Since there are 1/2 d openings per unit length, the leakage rate per unit 
area of the bladder if we assume a unit length for the depth is 


W' = — 


W_ 

2d 


pD c 


,mra..2 . , /nirhd 

2 d 3 V' Sln ( ~d~ ) tanh ( ~d~ } | 
h + 3 2 2-* 3 J 

•- it a n=l n 


(2.4.13) 
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A plot of W'/pDc as a function of a/d is shown in figure 2.8. 
v o 

As shown a/d is equal to the percent open area. For these calculations 
h was fixed at 5 rails and a at 1 mil and d was varied. This is 
unimportant, however, since W/pDc o depends only on the ratio a/d 
(as well as h/d, h and d) . From the plot we draw the following 
conclusions . 

o For small percent open area the leakage rate increases 
linearly with a/d and very rapidly. At 15% open area the 
barrier permits about 60% of the leakage rate without any barrier 
at all. Thus a barrier which has more than a few percent open 
area will not be effective. To be really effective a barrier 
must be about 99% intact. 

o For open areas greater than about 15% the approximately linear 
dependence of the leakage rate on slit width no longer holds. 

This is apparently the region in which adjoining slits begin 
to significantly influence each other since the leakage rate 
for one slit is proportional to the slit width for large 
widths. This agrees in order of magnitude with our previous 
results for circular openings where it was concluded that 
about 10 hole radii separation was sufficient for non- 
interference. 
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a/d - % Open Area 

Figure 2.8 leakage rate for Barrier with Periodic Slit Openings - A Specific Case 





2.5 Effect of Hole Shape on Leakage Rate 


13 


Pin holes as well as cracks which appear in the barrier of a 
laminated bladder structure will not, generally speaking, possess the 
geometrically regular shapes for which an analysis is possible. It 
is necessary, therefore, to determine as well as is possible the 
effect that hole shape itself will have on the leakage rate. To do 
this we compare the leakage rate through holes of different shape but 
having the same area. It seems reasonable to expect that the leakage 
rate will be about the same for holes having roughly similar shape 
but will diverge for holes which are very dissimilar , for example, a 
thin crack compared with a circular opening. The results obtained 
bear out this assertion. 

First a comparison of some results obtained using 
f 9) 

program DAP comparing circular and square holes are 

Table 2.1 

Hole Size Leakage Rate 

Side or Diameter Square Hole 

2 mils 1.0 x 10 -12 

8 mils 3.7 x 10 

It can be seen that the leakage rate through a square hole and a 
circular hole of equal area are close enough to each other so that for 
practical purposes they may be considered the same. The above calcula- 
tions were obtained for a bladder consisting of two layers of Teflon 
with Nitrogen gas diffusing. The results, of course, depend only on 
the geometry and not on the particular solute chosen. 

A further comparison can be made using the analytical results of 
sections (2.1) and (2.2) for circular holes and retangular slits. A 
value is chosen for the hole radius, a^ and the leakage rate obtained 
for typical conditions using the curves for circular holes. A slit 


computer 
given below 

gms 

2 

sec(mil) 
Circular Hole 
-12 

1.1 x 10 x 
-13 

3.9 x 10 
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is chosen of width a g and length £ which has the same area as the hole 

and the leakage rate through this structure obtained from the appropriate 

curves (Section 2,2). In order for the slit results to be applied to 

a rectangle it is necessary to choose £ >> a g . The results are shown 

in figure (2.9). The solid line gives the leakage rate through the 

circular hole and the symbols denote the leakage rate through rectangular 

slits having length to width ratios of 2, 5 and 10. The values for 

£ = 2a are somewhat low since these rectangles are not long enough to 
s 

be considered as infinite slits for which the analysis is valid. This 

is probably also true of the £ = 5a rectangles but less so. The leakage 
^ s 

rate for the £= 10a case should be accurate to within a few percent. 

s 

Figure (2.9) not only gives a comparison between rectangles of 

various shapes and circles of equal area but also between rectangles 

of equal area and different ratios of length to width. It appears from 

the figure that there is only a weak shape dependence for holes up to 

about 2 mils In radius and rectangles of up to about the same width. 

For larger holes the longer rectangles, £ = 10a , appear to be diverging 

s 

rapidly in leak rate from the value for circular holes. It appears, 
then that for holes (^ 1 mil radius) the shape could be changed from hole 
to square and to a long rectangle (£ = 10a g ) without changing the leakage 
rate appreciably so long as the area remains fixed. For larger holes 
there would be a greater variation in the leak rate. However, even for 
the 3 mil radius hole the change to a rectangle having £ = 10a g would 
only change the leak rate by about 20% (increase). 
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3.0 STEADY STATE LEAKAGE THROUGH HOLES IN METAL FOIL BARRIERS 

The leakage rate through a laminated bladder structure with a 
circular hole in the barrier was calculated in section 2.1. In this 
section we calculate the rate of leakage of pressurant gas into liquid 
propellant through a circular opening in an uncoated barrier, that is, 
a barrier consisting only of metal foil. The leakage rate through 
such a structure has intrinsic interest since this is the simplest 
type of bladder. The results also afford the opportunity to evaluate 
the effectiveness of the Teflon coating as an inhibitor of leakage. 

3.1 Leakage Through a Circular Hole in a Metal Foil 

Separating Pressurant Gas From Liquid Propellant 

A finite body of liquid propellant is separated from the pressurant 
gas volume by an impermeable barrier containing one small circular hole 
of radius a. The total pressure is the same throughout the pressurant 
gas-liquid propellant container so that there is no net differential 
pressure across the barrier. Assuming that the hole is formed after 
the pressurant and propellant are in place, after formation the 
pressurant gas will be dissolved in the liquid within the radius of the 
hole and will subsequently diffuse into the body of the liquid. 
Simultaneously propellant is vaporized into the pressurant gas there and 
subsequently the vapor diffuses into the body of the pressurant gas. 
Equilibrium occurs when the total body of liquid is uniformly saturated 
with dissolved pressurant gas and the pressurant gas everywhere contains 
propellant vapor with a partial pressure corresponding to the vapor 
pressure of the liquid at the system temperature. 

If the hole in the barrier is small and there are large volumes of 
gas and liquid present (relative to the dimensions of the hole) , it 
should be possible to obtain a reasonable estimate of the leakage rate 
by considering the system as infinite with mass flow rate by diffusion 
between two semi-infinite media connected by a small opening. It will 
be shown subsequently (Section 4.3) that the mass flow rate obtained 
in this manner is a good approximation to the actual leakage rate between 
two finite volumes for most of the time span during which mass flow takes 
place. 
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Gas and Vapor 


Liquid and Dissolved Gas 


Figure 3.1 shows the geometry of the problem to be considered. Con- 
sidering the vapor-pressurant gas mixture as a perfect gas the partial 
pressure of the pressurant gas in the gas phase (x > 0) obeys the equation 
(at uniform temperature) 


9p. 0 

— ~ = D V p, x>0 (3.1.1) 

3t g H 1 


Within the liquid, the dissolved pressurant gas will constitute only 
a small fraction of the mass of the mixture. Thus if ^ is the concen- 
tration of pressurant gas in the liquid propellant, c^ = Pg/( p £ +P g) 
p /p , the mass fraction c^ obeys the equation 


3c 2 

3T " V C 1 


X < 0 


(3.1-2) 


We have designated by the diffusion coefficient of pressurant 

gas through the gas-vapor mixture and denotes that for the diffusion 

of dissolved gas through liquid. The quantity p refers to the mass of 

T O 

dissolved gas per cm in the gas liquid mixture, p^ the mass of liquid 
per cm 3 there and p = p + p is the mixture density. 

The equations (3.1.1) and (3.1.2) must be solved subject to the 
following boundary conditions. There must be no leakage of gas across 
the impermeable boundary for r > a. Thus for r > a and x = 0, since 
the mass flux vectors are proportional to the gradients of p^ and c^. 
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(3.1.3) 


3c 1 

9x 9x 


For r < a (and x = 0) the mass flux must be continuous. In the liquid 

the mass flux vector is - pD Vc- and its x component is -pD (9c-/9x) i 
^ £ J. 1 

where i is a unit vector in direction x. In the gas phase the mass flux 
vector is -m^ D^Vn^,with x component, -m^ D^(3n^/9x) i . Here, m^ is 
the mass of the pressurant gas molecule and n^ is the number of such 
molecules per unit volume. The number of molecules per unit volume is 
related to the partial pressure (for a perfect gas mixture) by 


n l ” RT P 1 


when A is Avogadro's number, R is the gas constant in appropriate units 
and T is the absolute temperature. The mass flux in terms of the partial 
pressure is then 


A 8p l 

- m i ° g It aT 


i = 


M 1 D 
. 1 B 
RT 



(3.1.4) 


where = m^A is the molecular weight of the pressurant gas. Thus 
conservation of mass requires 


3c- M-D 3p- 

pD — — = -^-1 
p % 9x RT 3x 


(3.1.5) 


for x = 0, and all r since both 9p^/3x and 3c^/9x vanish for r > a. 

Finally for r < a and x = 0, we assume local thermodynamic 
equilibrium between the pressurant gas dissolved in the propellant and 
the pressurant gas dissolved in the gas phase. This gives the result, 
valid for dilute solutions of pressurant in propellant, 


c l = < p^ (3.1.6) 

when k depends only on the temperature. Here k is, of course, the 
Henry's Law Solubility for pressurant gas in liquid propellant. 
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In a cylindrical coordinate system, and for the steady state, 
eqs . (3.1.1) and (3.1.2) take the form 


, 1 3p l 
8r 2 r9r 


+ 



X > 0 (3.1*7) 




+ 



= 0 


x < 0 (3.1.8) 


A solution to these equations subject to the boundary conditions 
enumerated above, eqs. (3.1.5) and (3.1.6) must be found in order to 
calculate the steady state leak rate through the opening. The justifi- 
cation for obtaining a steady state solution is our surmise that the 
transient in the neighborhood of the opening should be short-lived so 
that the leakage rate will be virtually equal to the steady-state rate 
over most of the time span of interest. It will be shown later that 
the transient is indeed short compared to times of interest to a typical 
space mission. 

For large values of V r 2 + x 2 the pressure p x must go to the total 
system pressure, p Q , say, since there the gas will contain no vapor. 

So we put ip = p Q - Pj^ for x > a and 

1JL + I It + i-t = 0 (3.1.9) 

. 2 r 3r ~ 2 
dr 9x 

The function ip tends to zero for large values of r or x. Taking the 
Hankel Transform of order zero of eq. (3.1.9) with respect to the 
radial coordinate r, defined as 


iKk,x) 


f°° 

r t|j(r,x) J Q (kr)dr 
o 


(3.1.10) 


gives as the differential equation satisfied by (k,x) 

d2 *fe» x > - - k 2 *(k f x) = 0 
cx 


(3.1.11) 
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The solution to this equation which vanishes for large x is 


if)(k,x) « A 1 (k) e 


A^(k) being some function of k which must be found from the boundary 
conditions. From the inversion theorem for Hankel Transforms 


:,x) = f k A 1 

* n 


(k) e ^ J (kr) dk 
o 


(3.1.13) 


and therefore 


P 1 (r,x) = p Q - k A^C^e J o (kr) 
J o J 

x > 0 


(3.1.14) 


The solution for c^(r,x) is the same except for the sign of the 


exponent 


c^(r,x) = k A£(k) e J o (kr) dk x < 0 (3.1.15) 

■* o L 


Applying the condition of conservation of mass flux at x = 0, 
eq. (3.1.5) leads to 


E3 



M, D 

§m 

.1 8 

■ 

RT 

B 



which, since 

i 


B 



1 o 


A 2 (k) = RT ^ A l (k) 


Z o 


(3.1.17) 


The condition of local thermodynamic equilibrium, eq. (3.1.6) gives, 
after using eq. (3.1.17) 


K p - k A (k) J (kr) dk = — 

O 1 O ivl 

L J o 


M. D r°° 


— k A- (k) J (kr) dk (3.1.18) 
pD„ J 1 o 
K l J o 
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which must hold for r < a. Finally for r > a, 9p^/3x = 0 at x = 0, 
or 


k 2 A x (k) 


J (kr) dk = 0 
o 


(3.1.19) 


Hence we have the following set of dual integral equations for the 
unknown function A^(k). 


k A 2 (k) J Q (kr) dk jj“D 


+ I 

RT PD £ k 


(r < a) (3.1.20) 


k A (k) J (kr) dk - 0 
1 o 


The solution to this set of equations is known ^ ^ . 


A x (k) = r 


_ 2 / ^ sin ka \ 

" ir ( M D . , 2 J 

\l + _i _S_ I k / 

' 1 nm _ tv 


RT pD £ k 


r > a (3.1.21) 


(3.1.22) 


Combining eqs. (3.1.22), (3.1.17) and (3.1.15) gives for the concentra- 
tion of dissolved pressurant gas in liquid propellant, 


c x( r > x ) = T 


_ 1 V / P o \ 

'" RToD A 1 + ^vJ 


. g . 1 P ..M e kx J Q (kr) dk (3.1.23) 


RT pD £ k 


From this result the leakage rate per unit area is 
w(k) = pD. 


9C 2 

= 2 

M i 

/ P o 


9x 

7T 

RT °g 

1 M i 

D ) J 


x=o 

'1 + — 
RT 

_a !/ J 
pD * K 


sin ka J (kr) dk 
o 


(3.1.24) 


The integral vanishes for r > a as it should. Its value for r < a is 

. ,Tl 2 

l//a - r . 


Thus 


w(r) 


2 M l D g / P o \ L 

' ET + 


(3.1.25) 


RT pD^ k 
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The total leakage rate is 


W = 2tt w(r)r dr 
o 


(3.1.26) 


and is given by 


w ■ 4a ir ( + M^V i ) 


(3.1.27) 


RT pD £ K 


The leakage rate thus depends directly on the total system pressure 
p^, the radius of the hole, a, and the group of parameters D M^/RT. 

The retarding effect of the solubility in and diffusion of the gas 
through the liquid is contained in the group of parameters (M^/RT) 

(D /pD ) (1 /k) , which can be recognized as the ratio of the gas to the 

liquid permeation rates. We estimate the value of this ratio for the 
typical case of Helium diffusing through ^0^. For Helium at a tempera- 
ture of 273°K, M /RT - 1.78 x 10 According to the data of Chang 
* _ ^ _•£ 

and Gokcen(12) , k ^4.3 x 10 atm , The density of N O is about 

3 ^ ^ 

1.45 gms/cm . No data is available for the diffusivity of Helium in 
liquid ^0^, but the value should not exceed 10 ^ cm^/sec^. For 
gases through gases, the diffusivity is always of order 1.0. Using these 
numbers we find 

iVl in 6 

RT pD^ k * 10 

We conclude that the leakage rate is controlled by the solubility and 
diffusivity of the gas through the liquid and for all practical purposes 
eq. (3,1.27) reduces to 


W = 4a p k p Q D £ (3.1.28) 

Equation (3.1.28) is the same result that one would get if the 
gas phase was entirely ignored and the problem was considered to be 
leakage through a hole into a half-space containing the liquid propellant, 
the pressurant gas within the radius of the hole being held at the value 
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According to eq. (3, 1.28) the leakage of pressurant gas into the 
liquid is directly proportional to the solubility of the gas in the 
liquid, the total system pressure and the diffusivity of the gas in the 
liquid as well as the 'radius of the hole. In the next section the 
leakage rate through such an opening will be compared with that obtained 
for leakage through a similar hole in the barrier of a Teflon laminated 
structure . 


3.2 Effect of the Teflon Coating on Leakage 


Assume that the hole in the barrier Is 1 mil in radius and is 

circular. Take the total system pressure as 40 atmospheres. The value 

—5 2 

of D is estimated at 10 cm /sec. According to Jost(l) all diffusion 
& 

coefficients in liquids have a value in the neighborhood of this number. 
In the case of system gas diffusing through ^0^, the solubility, 
k, is somewhat larger than for the system He-N^O. . For N- gas In ^0^ , 

K = 1.8 X 10 atms (1). The leakage rate obtained for the two systems 
from eq. (3.1.28) is shown below 


System 

He-N o 0, 
2 4 


H 2 - N 2°4 


W gms/sec 
2.5 x 10 -11 


1.06 x 10 


-9 


For purposes of comparison the leakage rate through a one mil radius 
hole in the barrier separating 5 mils of FEP on the gas side from 5 mils 
of TFE on the liquid side at the same total system pressure Is shown 
below. 

Pressurant Gas W gms/sec 

He 3.0 x 10 -11 

N 2 1.67 x 10' 11 

Thus the Teflon coating reduces the leakage rate by about two orders 

of magnitude in the case of N 2 but does not significantly affect the 

leakage rate in He-N 2 0^ systems. The reason for this is that the diffu— 

sivity of He into Teflon is the same order of magnitude of He through 

N„0, . These conclusions as stated above are based on estimated values of 
2 4 

the diffusion coefficients of both gases [He and N 2 through liquid N^] 

c 2 

of about 10 cm /sec and therefore must be regarded as tentative since 
at present there is no experimental data for the diffusivity of N 2 and 


H g gases into liquid N 2 0^. 
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4.0 TRANSIENT LEAKAGE THROUGH BLADDER STRUCTURES 


m 




m 


It was pointed out in section 1.1 that for a closed system, that 
is a gas reservoir separated by a permeable bladder from a liquid 
reservoir, the only steady state solution to the diffusion equation is 
the trivial one corresponding to an everywhere uniform concentration 
and no leakage. What must occur in reality is that after loading the 
system with pressurant and propellant, or, if after loading, a rupture 
In the diffusion barrier occurs there is a transfer of mass between 
the two compartments until the liquid is everywhere saturated with gas 
and conversely. This process does not take place at steady state but 
is accompanied by a constantly diminishing leak rate as each tank nears 
saturation since the driving force for the mass transfer steadily 
decreases. 

However, if the leakage rate through the bladder is slow, that is, 
the holes in the barrier are small and if the pressurant and propellant 
reservoirs are relatively large, that is, so large that at the given 
leak rate the time required for saturation is long, it is possible to 
consider quasi-steady state solutions to the leakage problem. The steady 
state leakage rates calculated in sections (2.0) and (3.0) apply to this 
situation. This type of solution will give an accurate estimate of the 
leakage rate for periods of time over which the concentration of the 
solute at the bladder surface does not change appreciably. If the 
calculations are to be applicable to a real bladder situation it is 
clearly necessary that the transient time of the bladder structure 
itself be short compared to the time span over which the solute concen- 
tration in either the liquid or gas reservoirs changes appreciably. The 
latter time span must be large, of course, if the bladder is functioning 
properly, that is, preventing leakage. One of the main purposes of the 
transient calculations, therefore, is to determine if the relaxation time 
for the bladder is short compared to that for the system as a whole. 

If the ratio of the two times is sufficiently small the leakage rate 
through a bladder structure can be calculated from the steady state 
results since the integrated leakage over the bladder transient will be 
small compared to the steady state leakage integrated over periods of 
time of the order of a mission lifetime. 
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4.1 Approximate Calculation of the Transient Period For a Laminated 

Bladder with a Circular Hole in the Barrier 

Consider a semi-infinite slab of material 1, extending from x = 0 
to x = h axially and to infinity radially. The plane x = 0 is covered 
by an impermeable sheet except for an opening of radius a which is 
exposed to the diffusing vapor. The concentration on the plane x = h 
is held at zero. Initially, the concentration of the diffusing sub- 
stance throughout the slab is zero. A cross section of the bladder 
is shown in the schematic. 



The basic difficulty with all problems of this type arises due to the 
mixed boundary condition which must be satisfied on the plane x = 0. 
Thus, one specifies the concentration for r < a and the gradient of 
the concentration (zero) for r > a. If it were possible to specify 
either the concentration or the gradient over the whole plane x = 0 
then the problem would be straightforward. This cannot be done exactly, 
of course. As an approximation, however, one can postulate that the 
flow rate across the opening (r < a, x = 0) experiences a very short 
transient before it takes on the final steady state value. The 
difference between the leakage rate, which is the flow across the 
plane x = h at any given time and its steady state value is then due 
to the storage of diffusing material in material 1. This approach 
might be expected to give accurate results after the passage of a 
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period of time short compared to the relaxation time for the bladder as 
a whole because equilibrium should be reached very rapidly in the 
immediate vicinity of the hole. 


To find the transient behavior of the leakage rate under these 
assumptions we solve 


2 ^1 1^1 A l_Ifl 

v C 1 2 + r 3r + 2 D 3t 

dr 8x 1 


(4.1.1) 


Subject to the boundary and initial conditions. 


c^r.h.t) - 0 
c 1 (r,x,0) = 0 


-D 1 g^-(r,0,t) = f(r,t) 


where f(r,t) is zero for r > a, and is equal to the steady state flow 

2 

rate divided by it a for r < a. That is 


f (r , t) = — 2 


r < a 


f (r, t) = 0 


r > a 


The Laplace Transform of Equation (4.1.1) then gives, using (4.1.2b) 


2 2 
3 c, (r,x,s) . 3c.(r,x,s) 3 c..(r,x,s) 

t c i <r ' I '- s) ■ 2 — + r — + a — 

1 dr 3 x 


where 


c.(r,x,s) 


e St c^(r,x,t) dt 


Taking the zero order Hankel Transform of Equation (4.1.3) gives 


where 


d c. (k,x,s) , 

z (k + ~) c- (k,x,s) = 0 

dx 1 


r°° 

c 1 (k,x,s) *> r J (kr)c (r,x,s)dr 


(4.1.4) 
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r. i 


: tl 


The solution to equation (4.1.4) which satisfies the condition (4.1. 
is 


c. (k,x,s) = A (k,s) 


sinhVk^+-^- (h-x) 
1 


sinhv / k2+~- h 


(4.1. 


; m 


The Hanke 1-Lap lace Transform of the conditions (4.1. 2d) and (4.1.2e) 
determines A^(k,s) , that is 


dc 1 (k,0,s) ra 

-D — 51 J (kr> dr j 2 

J o J o Tta 


r -st 
L e 


dt 


W , aJ, (ka) 

O 1 1 


Tra 


2 s 


(4. 


j g§ 

I L. i 


We thus find 


c 1 (k,x,s) 


... sinh k 2 +j^- (h-x) 

W a J, (ka) , D, 

O 1 ± 

2D, k s 
tra 1 


V? 


s 

D, 


+ ~ cosh\k2+^- h 


(4. 


The Hankel-Laplace transform of the flow rate through the plane x 


W(k,x) * -D 


dc^Ckjh^s) W q a J^(ka) ^ ^ 


1 dx 


?ra 


2 k 


coshs k^+j^- h 


(4. 


Performing the Hankel Inversion gives 


W 


W(r,s) =-ff 


r°° J (kr) J, (ka) 


tra 


o cosh'vk^l* — h 
D 1 


dk 


(4. 


1 

The total leakage rate, or rather its Laplace transform is given by 


/**> 2W f<>» r°° J (kr)J 1 (ka)dk 

W(s) - 2tr r dr W(r,s) = — | f j r dr — _ — — 

•'o a " ■'o ■'o cosh v k2+ — — h 


(4. 


2a) 

5) 

1 . 6 ) 

1.7) 

h is 

1 . 8 ) 

1.9) 

1 . 10 ) 
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Interchanging the order of integration, and using the delta function 
character of the integral with respect to r, i.e., 

f r V»> d. - ^ 

; o 

where Sflc) denotes the Dirac delta function, we get 


W(s) 



(4.1.11) 


Performing the Laplace inversion, 
yields 



which is straightforward in this case 


-(2n+l) 2 

e 



(4.1.12) 


The leakage rate as a function of D^ has been plotted in Figure (4.2). 
From the figure, it can be estimated, that for a 5 mil thick slab that 

the transient period will last about 700 seconds if D is as large as 

fi i —9 X 

10 cm /sec, while for T> 1 as small as 10 this period will last 194 

hours . 


The value of the diffusion coefficient for pressurant gases through 
—7 2 

Teflon are about 10 cm /sec (see the experimental section of this 
report) . From figure (4.2) the transient time for a 5 mil thick Teflon 
laminated structure is about 3 hours or 10 seconds. For a bladder 
composed of two Teflon Laminates the time required would be about 4 times 
as long since the transient time depends on the square of the bladder 
thickness (eq. 4.1.12). 

Although the above analysis is based on an approximate model, it 
is expected to give the accurate asymptotic solution because for long 
times the leakage rate across the opening must be the steady state 
leakage rate. This has been verified by running the problem on the 
digital computer using DAP. Figure (4.3) shows the analytical results 
for a 5 mil thick Teflon slab and a barrier with a one mil radius 
circular hole and the results obtained using the digital computer. 

It can be seen that the results are in fair agreement even over the 
early part of the transient and converge nicely for long times. 


4-5 


cm Teflon thickness 



4-6 


Figure 4.2 



mm 



seconds 











fra 

w 


m 


The following general conclusions can be drawn from the above 
results . 

o Transient times for a typical laminated bladder structure are 
of the order of a few hours. This is very short compared to 
the mission lifetime required for deep space exploration, 

o For a one mil radius hole in the barrier the steady state leakage 

rate for N 9 gas through two layers of Teflon each 5 mils thick 
^ -12 

(see section 3.2) is about 3,42 x 10 gms/sec. In a period 

of 10 5 seconds which should greatly exceed the transient time 

of the bladder structure, the accumulated steady state leakage 

rate would be 3.42 x 10 - ^ gms. The accumulated transient 

leakage rate would be less than this as can be seen from 

Figure 4.2. Such a low integrated leakage would scarcely 

affect the concentration of N ? gas in the liquid propellant 

z 12 

over a period of time of the order of 10 seconds which is 
very long compared to the 10^ second or so bladder transient 
time. 


* 5 ? 
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4.2 Further Remarks on the Transient Problem 


Another instructive example for which an asymptotic (large values 
of the time) solution can be obtained is that of the transient diffusion 
of pressurant gas through a circular hole in a metallic foil bladder 
into an infinite body of liquid. This solution is developed in detail 
in reference 13. Here for the sake of brevity only the results will be 
given. 


The initial conditions on the problem are that there is no gas 
dissolved in the liquid at time t * 0, and it is assumed that the concen- 
tration of the pressurant gas within the opening is held at the constant 
value c q . This latter assumption is justified in section (3.1) and is 
due to the fact that diffusion through the gas phase is much more rapid 
than through the liquid. The leakage rate in gms/sec through a hole of 
radius a is then found to be 


W(t) = 4apD„c + a 
t o 


3 p Vo 

V 


. 33 5 pD £ C o 

+ 48 a 2 

<v> 


(4.2.1) 


This formula is valid for sufficiently large values of the time and 
clearly goes to the steady state leakage rate 


W « 4apD.c 
o £ o 


(4.2.2) 


for large values of the time. 

From eq. (4.2.1) the leading term in the time is small compared 
with the first when 




« i 


or when 


4D„t » a 
£ 


(4.2.3) 


For a one mil hole, and an assumed value of of 10 cm /sec, the 

leading transient term in eq. (4.2.1) is 1% of the steady state term when 

t - 16 seconds. At this time the ratio of the next transient term to the 

-4 

steady state term is 1.35 x 10 , i.e., very small. Thus the transient 

is very short-lived for the hole sizes under discussion here. 


U4 


Let us now compare this transient time with the time required to 
saturate a given quantity of propellant by leakage through a one mil 
radius hole at the steady state rate, eq. (4.2,2), At 40 atms the 
equilibrium concentration of He gas in N^O, is about 

/ O ^ H 

1.74 x 10 (gms He)/(cm J N 9 0,) according to the data of Chang and Gokcen. 
One ft of propellant contains 2825 cm and at saturation 0.492 gms of 

Helium. The steady state leakage rate from eq. (4.2.2) is 

—11 3 

2.5 x 10 gms /sec. The time required to saturate 1 ft of propellant 

at the steady state leak rate would then be 1.96 x lO 1 ^ secs or 

5.46 x 10^ hrs. Actually it would require longer since the steady 

state leak rate above is greater than the actual leak towards the end 

of the time period required for saturation. Thus the transient period 

is extremely short compared to the time required for the transfer of 

any significant quantity of pressurant gas. Furthermore the leakage 

rate can be accurately calculated from eq. (4.2.2) for periods of time 

of the order of the time required for a deep space mission. 

Finally it should be pointed out that the DAP program is now 
available for solving many transient problems involving two laminates 
of permeable material separated by a barrier with a hole in it. The 
hole shape may be circular, square, or rectangular. Running time for 
a typical transient problem, say leakage through two 5 mil layers of 
Teflon with a 1 mil circular hole in the barrier, will be from 10 
to 20 minutes. The program is not particularly well suited for problems 
like the one discussed above where the transient time is very short, 
however. The computing time required is too long. For more detailed i 
information on the program see reference 9. 


rj 
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5.0 INSPECTION TECHNIQUES FOR DETERMINING THE PRESENCE 

OF MICROSCOPIC HOLES IN METAL FOIL 

Before the leakage rate through a given laminated bladder structure 
can be predicted, the size, shape, and distribution of holes In the 
barrier foil must be known to sufficient accuracy* Of these quantities 
the size and distribution are most Important since as was shown in 
section (2.5) the leakage rate Is only weakly dependent on the hole 
shape for not too gross differences in the shape. For small holes, 
however, the leakage rate is strongly dependent on the size of the hole 
In the sense that the leakage per unit open area Is not constant. Thus 
the leakage rate through two circular holes will exceed that through one 
circular hole having an area equal to the sum of the areas of the smaller 
two holes. In addition, a certain minimum hole separation is necessary 
in order to be able to treat the holes Individually. According to the 
results of section (2.3) this should be about 5 times the major open 
dimension of the hole. Thus an accurate estimate of the leakage rate 
through a bladder structure can only be made if a considerable amount 
of detailed information is available describing the holes in the barrier. 

One of the most economical and well developed techniques for 
obtaining such information appears to be dye testing. This technique 
utilizes liquid dyes which are extremely good surface wetters along 
with subsidiary techniques for observing the hole. It appears that such 
techniques could be used for locating and photographing holes as small 
as 1 mil diameter in metal foil. 

One such process Is known as dy-chek and is a product of Turco 
Corporation. This test is carried out by first cleaning both surfaces 
of the foil carefully and then applying liquid dye-penetrant which is 
non-corrosive. After the dye— penetrant has dried an emulsifier is added 
followed by a developer. The presence of a hole is indicated by the 
appearance of dye coloring on the opposite side of the foil. An elabora- 
tion of this process is Zyglo, a product of the Magnaflux Corporation. 

It is similar to dy-chek except that the dye is flourescent and can be 
viewed and photographed under black light. It should be possible to 
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obtain both the shape and size of the openings from enlargements of such 
photographs. After measurement of the holes in this manner the method 
could be checked by a gas leak rate test which would give the total open 
area available for gas leakage. 

More exotic techniques can be imagined such as 3 particle attenua- 
tion, or optical methods using laser light. However these methods should 
not be attempted until the simpler and more economical methods suggested 
above have been exhausted. 
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6.0 SUMMARY OF ANALYTICAL RESULTS 


In this section we summarize those general conclusions which can be 
drawn from the analytical study, compare the analytically predicted leak 
rate with experimentally measured values and point out those areas which 
need further attention. 

6.1 Summary of Results 

As a result of the work reported In sections one through five above, 
and elsewhere, several general conclusions concerning leakage and its 
calculation through laminated bladder structures consisting of two layers 
of permeable material separating a barrier with a hole can be stated. 

The most important of these for an understanding of the overall process 
are 

o The transient time for a typical bladder structure is small 
compared to the time required for a significant change in the 
concentration of the gas In the liquid propellant. This means 
that meaningful leakage rate calculations can be made by assuming 
the bladder structure is exposed to steady state conditions on 
each of its surfaces. 

o Concentration gradients in the gas phase can be ignored. This 
is shown in Appendix A. The reason is the gas offers essentially 
zero resistance to diffusion compared to the Teflon. 

o Concentration gradients in the liquid phase (propellant) can 
also be ignored to a first approximation. This conclusion, 
however, must remain tentative at present since there are no 
experimentally measured diffusion parameters available for the 
pressurant gases (N 2> he) in the common earth storable propel- 
lants (n 2 o 4 , n 2 h 4 ). 

o As a consequence of the above results the upper bound for the 
leakage rate through a bladder structure can be calculated by 
specifying a constant concentration on the solute side, zero 
on the opposite side of the bladder. The solute here designates 
the diffusing component. 
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o The shape of the hole in the bladder only weakly affects the 
leakage rate per unit area through the hole over a wide class 
of hole shapes. For example, a rectangular slot having a length 
10 times its width will permit about the same leakage rate as 
a circular hole of the same area. 

o Openings in the barrier can be treated as independent provided 
they are separated by about 10 or more times the major dimension 
of the opening. 

o Provided the diffusion coefficients of the pressurant gases 
(N^, He) in liquid propellant (N^O^,, ^H^) are no smaller than 
about 10 ^ cm^/sec, the Teflon coating on the barrier foil 
reduces the leakage rate through the bladder structure signifi- 
cantly over what it would be for an uncoated bladder (plain foil 
with a hole) . 

6 . 2 Comparison of Ca lculate d with the Experimentally Measured Leak Rate 

i 

To determine the validity of the analytical approach experiments 
were carried out to measure the leakage rate through a typical laminated 
bladder structure containing circular holes in the barrier* The sample 
bladder structure consists of a 10 mil layer of FEP separated from a 10 mil 
layer of TFE by a one mil thick aluminum foil. The foil was drilled with 
100, 5 mil radius circular holes in a square pattern, approximately 
125 mils on centers. The sample laminated bladder structure was obtained 
from the Dielectrix Corporation; the holes were previously drilled in 
the foil by National Jet Corporation. 

The tests were carried out by initially evacuating both sides of the 
test assembly, then pressurizing one side with Helium gas, the other with 
N^O^ vapor. The two sides are separated by the laminated bladder struc- 
ture. To eliminate the leakage during the transient period from the 
measured value of the leak rate, the gas and vapor were flowed through 
their respective sides of the assembly for several hours before the system 
was closed off. After close-off sufficient time was allowed to elapse 
to obtain a measurable quantity of the solute in the opposite chamber. 

The samples were analyzed using gas chromatography. 
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The experimental runs were carried out with the results shown in the 
table below. The leakage rate through the assembly was ca^ 

Ca.v^7 'if , . , 

the method oUt4ijtad" i r ^ the, .d es ign- guid e using the - 

qh rveo - for 'leakage th r ou g h - cifycul gfc * 



la atsi 

enerated 
/ry^_ 


Test 


1 


2 


3 


Temp . 

Press. 

(Atm) 

Duration 

(hrs) 

W gms/sec 
(Exper.) 

W gms/sec 
(Calc) 

187 

9.6 

2 

He-*N„0. 2.72 x 10~ 9 
2 4 

4.20 x 10 -9 




N„0.->He (Not Measurable) 
2 4 

200 

10 

64.5 

He^N_0. 15.5 x 10 9 
2 4 

4.8 x 10 -9 




N o 0.^He 3.1 x 10~ 8 
2 4 

4.35 x 10 -8 

150 

4.9 

44 

He-^N-0. 17.35 x 10 -9 
2 4 

2.11 x 10 -9 




N„0. ->He 1.6 x 10 -8 
2 4 

1.7 x 10 -8 


All of the tests were carried out with the TFE on the ^0^ side. 

It was subsequently found that Test 1 was invalid because the gas 

and vapor were not allowed to flow through the system for a sufficient 

length of time prior to closing the system off. The measured leakage 

rate is therefore an average value over a transient period and for this 

reason is low. The test was also of too short duration to allow the 

determination of the N_0. rate- 

2 4 

For tests 2 and 3 the system was not only exposed to Helium and 
for a sufficient length of time before close off but also the test 
was allowed to run for a period long compared to the estimated transient 
time for the laminated bladder structure. In this way valid determina- 
tions of the steady state leak rate were obtained. 


The calculated values u f the leak Tale were ublallied using the * 
sel ubillty and dlffuSlvity dat - a of tabloc 9-1, 9 - 4 t —and— 9~-7 . -The— 
d rif fu sdrvit y dat a used was that obtained from-uuldirect i o nal ex p eriment s-. 
The experimental leakage rate was obtained by measuring the total accumu- 
lated leakage and dividing by the elapsed time. From the results of 
tests 2 and 3 above, it can be seen that the predicted Helium leak rate 
is consistently lower than that measured while the predicted ^0^ rate 


Jr- /J/hS/9 ' ' #*4~*** 7 
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is reasonably close to the measured value* However, there are several 
factors which must be considered before evaluating the accuracy of the 
predicted leak rate* The most important of these are the following. 

(1) After completing the tests the laminated bladder structure was 
removed from the test assembly for examination. It was found that the 
TFE which was on the ^0^ side of the bladder was completely separated 
from the aluminum foil. The FEP bond remained intact however. In 
addition, there were cracks in the foil barrier at various places along 
the surface where the bladder met and was compressed by the flanges of 
the test assembly. Both these defects tend to increase the measured flow 
rate over what it would be for an unfaulted bladder structure. The 
theoretical calculations were based on the assumption that the bladder 
material was not separated from the foil barrier and that the only holes 
in the barrier were the drilled circular holes. Quantitatively the 
effect of the separation and rupture of the barrier on the leakage rate 
is difficult to account for precisely since the open area of the crack 

is not known, and the amount by which separation increases the leak rate • 
is difficult to determine. Separation of one layer should not increase 
the leakage rate by more than 50% however, and a visual examination of 
the cracks Indicate that the total open area was considerably less than 
that of the 100 drilled 5 mil radius holes. Thus one would expect the 
measured leakage rate through a bladder structure which did not experi- 
ence separation and which had no cracks to be no less than about half 
the measured values reported above. In short the observed failures of 
the experimental sample could cause only a small change in the leakage 
rate, not order of magnitude changes. 

(2) The permeation rate for Helium through Teflon in the presence of 

counterflowing is higher than that for unidirectional flow of Helium 

alone. This result was observed experimentally by comparing the results 
of unidirectional and counterflow tests .and the results are-^. 
rtec o ide d When these values of the permeation rate are 

used, the following calculated leakage rates are obtained for tests 2 
and 3. 
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Test 


W(gms/sec) 

2 12.0 x 10 -9 

3 11.0 x 10 -9 

When these results are compared with the experimental leak rates given 
above, it can be seen that the experimental (as measured) and calculated 
leak rates are in fair agreement. When allowance is made for the 
observed defects in the bladder as discussed in (1) above, it appears 
that the calculated and experimental gas leak rates are in excellent 
agreement. Unfortunately, it is not possible to state this without 
reservation since no quantitative measure of the effect of the defects 
in the bladder structure is available, 

(3) The experimentally determined effect of counterflow on the 
leakage rate is to reduce it by about a factor of 10, (Se e — figure 

— whir l! sumina rirggs^the ex - pc riment a^ r ^rniTTrHh Thus when counter diffusion 
of Helium is taken into account the calculated leakage rates shown in 
the table for tests 2 and 3 would be reduced by about a factor of 10. 

Thus the agreement would not be nearly as good as is indicated by the 
table. On the other hand the experimental leakage rates are too high 
for the reasons discussed in (1) but not by as much as a factor of 10. 

(4) The experimental data is itself uncertain perhaps to as much as a 
factor of 10. Further, there is considerable scatter between different 
samples of Teflon. Since the experimental data used to make the theoreti- 
cal calculations was obtained from different Teflon from that used for the 
experimental leakage tests, there is some uncertainty involved. 

When all the possible sources of error such as those discussed in 
(l)-(4) are taken into consideration, it seems reasonable to assert that 
the prediction of the leakage rates to within an order of magnitude is 
acceptable. Our experiments and analysis indicate that the unidirectional 
diffusion data is adequate for this purpose and that the leakage rates 
calculated using this data is in satisfactory agreement with the experi- 
mentally measured leakage rates, 



6,3 Suggestions for Further Study 


The results obtained under this contract and reported above are all 
obtained by isolating the bladder structure from the rest of the system 
and calculating the leak rate as though it were in a fixed environment 
with respect to the diffusing components. It has been shown that so far 
as the gas space is concerned this is a. very good assumption over long 
periods of time (times of the order of that required for deep space 
missions). However, the assumption of a vanishing concentration of the 
diffusing pressurant gas on the face of the bladder exposed to the liquid 
propellant is subject to question. True, if the ratio of the permeability 
of the gas through the liquid is much greater than that through the Teflon 
the mathematical problem can be treated as though the concentration in the 
Teflon vanishes at the surface exposed to the liquid (Part 1, Appendix A). 
We do not know that this condition is adequately satisfied since we do not 
know the permeation rate of pressurant gases through liquid propellants. 
Even if it is, the model depends on the bladder surface remaining in con- 
tact with the liquid, that is, no vapor phase is formed in the propellant 

space. But according to Raoults* law (14) the pressure in the liquid 

space will increase due to absorption of the pressurant gas in the liquid 
tending to force the bladder away from the liquid and forming a body of 
vapor in the liquid space. In fact, the formation of the vapor phase on 
the liquid side of the bladder has been observed. 

The rate of formation of the vapor phase (we refer to the mixture of 

pressurant gas and propellant vapor as vapor phase) on the liquid side as 

a function of the degree of perforation of the diffusion barrier is a 
problem of considerable importance, since the presence of gas on the liquid 
side of the bladder can lead to engine malfunction if large quantities are 
present. Thus, the conditions under which a bubble forms and the rate of 
formation of the vapor phase for a given set of conditions should be 
studied. This can be carried out by extending the present study to in- 
clude the liquid propellant and the pressurant gas adjacent to the bladder 
structure and making allowance in the calculations for the motion of the 
bladder. Further experimental work to determine the diffusion parameters 
of the pressurant gases through liquid propellant will be necessary in 
order to carry out this study. 
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7 . 0 EXPERIMENTAL EFFORT 


The purpose of the experimental effort (Part II of the Bladder 
Permeation Program) was to support the development of an analytical 
model of bladder leakage by providing transport property data and 
identifying the mechanisms involved. In order to achieve the objective 
an experimental program was conducted to: 

1. Measure the diffusion coefficient and solubility of gases 
and propellant vapors in bladder materials. 

2. Assess the effects of two component counter diffusion. 

The experimental program consisted of tests with a 10 mil sample 
of polytetraf luoroethylene (TFE) , a 10 mil sample of tetrafluoroethylene- 
hexaf luoropropylene copolymer (FEP) and a composite sample consisting of 
1 mil of aluminum foil with 10 mils of. TFE bonded to one face and 10 mils 
of FEP bonded to the other face. The gases and vapors of interest to 
the program were helium, nitrogen, nitrogen tetroxide and hydrazine. The 
teflon specimens were purchased from Dilectrix Corporation and were 
fabricated by the spray method normally used for fabricating bladders. 

The propellants were from lots analyzed to meet MIL-specifications and 
the pressurant gases were industrial cylinder gases. 

For simple geometries, the permeation constant is directly propor- 
tional to the leakage rate. The mathematics involved and the mechanism 
are widely reported and discussed, for example. References 15 and 16. 
However, for the analysis of bladder leakage which involves flow through 
holes in an otherwise impermeable barrier, the transport property data 
required are solubility data in terms of a Henry's law constant and the 
diffusion coefficient. These data are seldom reported in the literature 
and, in fact, an extensive literature survey was reported in Quarterly 
Progress Report No. 2* which located only permeation constants. 

^Report No. 07282-6006-R000 dated 4 January 1967. 
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The articles included in the survey pertinent to the present study are 
References 17, 18 and 19 which give helium and nitrogen permeation con- 
stants in TEP and FEP. A continuing search for data, particularly for 
nitrogen tetroxide and hydrazine, yielded only four additional reports 
(References 20 to 23) all of which report permeation data for nitrogen 
tetroxide. Diffusion coefficients were not located for any of the 
systems of interest and only one source gave data which involves solu- 
bility in any of the systems of interest. Reference 23 reported that at 
75 +5°F after 183 days immersion in liquid nitrogen tetroxide, TFE 

increased 3.8% in volume and 2.6% in weight. This gives a solubility 

3 3 

of 56 mg of nitrogen tetroxide per cm of TFE compared to 22 mg/cm 
which could be predicted from extrapolating the low pressure tests 
discussed in Section 8.0. The comparison of the available permeation 
data with results of this study is given in Section 9.4. 

From an examination of the apparatus and results reported in the 
literature two types of experiments were selected to provide the required 
data. The first was an adsorption test to provide solubility data and 
the second was a permeation test to provide permeabilities and diffusion 
coefficients. These tests are discussed in Section 8.0 and the results 
obtained from the test program are given and discussed in Section 9.0. 
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8.0 TECHNICAL DISCUSSION 

The experimental program provided the measured values of the Henry s 
Law solubility constant, the diffusion coefficient and the permeability 
constant through the reduction of data from two types of experiments. 

The solubility of gases or vapors in a solid of known volume was measured 
by introducing a known quantity of each in a closed bomb and allowing 
time for equilibrium to be reached. From the equilibrium values the 
solubility was calculated. In those tests where the transient decay 
was sufficiently long, the diffusion coefficient also was calculated. 

The permeability constant for unidirectional gaseous diffusion through 
a membrane was measured by an apparatus in which a membrane was sub- 
jected initially to a small differential pressure. The differential 
pressure decreases in the closed system as a function of time due to 
permeation of gas from the high pressure to the low pressure side. From 
an evaluation of the differential pressure— time data, the permeability 
constant was calculated. 

The following sections describe the tests which were conducted as 
a part of this study. Section 8.1 describes the experiments using the 
adsorption apparatus and Section 8.2 describes the unidirectional 
permeation and counter flow tests using the permeation cell. 

8.1 Adsorption Tests 

In principle, the adsorption tests for gas solubility in solids 
is a relatively uncomplicated test that provides a maximum of information. 
In actual practice, it i9 capable of producing good results but requires 
careful calibration and attention to detail. The schematic diagram of 
the adsorption apparatus used in this program is given in Figure 8.1. 

The supply tank fill system, the evacuation connections and the associ- 
ated valving were omitted in order to provide a clear understanding of 
the basic apparatus. 

There are three key volumes which must be known accurately. These 
are: 1) the volume Included between the two valves which includes the 

transducer, 2) the volume of the empty adsorption bomb and all fittings 
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Pressure Transducer 
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Figure 8.1 Adsorption Apparatus 



and connections, and 3) the volume of the Teflon sample. In addition, 
the complete apparatus must be carefully temperature controlled for, in 
effect, it functions as a gas thermometer. 

The procedure which was followed during a test is the following: 

1. The sample volume and transducer section were swept with the 
gas or vapor and evacuated for a length of time sufficient 

to assure desorption of all gases. (Typically one to sixteen 
hours . ) 

2. A sample of gas or vapor was introduced into the calibrated 
volume between the two valves and its pressure was allowed 
to stabilize. 

3. The valve between the transducer and the sample was opened 
and the pressure in the system was periodically read on a 
potentiometer until no further decrease was observed. 

4. This valve was closed and a second sample of gas or vapor was 
admitted to the calibrated volume as in step 2. 

5. This second sample was added to the first as in step 3 and the 
procedure of steps 2 and 3 continued until the desired total 
pressure was attained. (Usually a total of four or five times.) 

The data were reduced in a direct manner when a permanent gas 
(nitrogen or helium) was used but by a more complicated technique when 
nitrogen tetroxide was used. For permanent gases, the mass (or 
standard volume) of gas in the system at any time was calculated from 
the pressure volume relationships. The difference between the quantity 
added and the quantity indicated at any time is the quantity adsorbed in 
the Teflon sample. 


CC(STP) gas 
CC Teflon 


273.15 
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( 8 . 1 . 1 ) 


where 
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Tef 
V B = 


TR 

^TR 

P 


EQ 


o Tr 

temperature, K 

3 

volume of the teflon sample, cm 

3 

volume of ads6rption bomb, cm 

3 

volume of transducer section, cm 

pressure in transducer section prior to addition to the bomb, atm 
equilibrium pressure in the bomb prior to the addition, atm 
equilibrium pressure in the bomb after the addition, atm 


When the nitrogen tetroxide adsorption data was reduced it was 
necessary to account for the dissociation reaction 


N 2 ° 4 ^ 2 N0 2 


( 8 . 1 . 2 ) 


which is a function of temperature and pressure. A second factor related 
to the fact that it was a condensable vapor. Liquid nitrogen tetroxide 
was held in the supply tank at the same temperature as the complete 
apparatus. When a vapor sample was drawn into the transducer section 
a variable quantity of liquid was entrained with the vapor. Thus, it 
was not possible from the initial transducer pressure-volume relationship 
to determine the quantity of nitrogen tetroxide added to the sample bomb. 
Since the adsorption of nitrogen tetroxide in Teflon is a slow process, 
it was possible to take a large number of pressure-time readings after 
the nitrogen tetroxide was introduced but before equilibrium was attained. 
From these measurements, it was possible to derive both the zero time 
pressure of the bomb (actually the gas density was used) and the diffusion 
coefficient of nitrogen tetroxide into the Teflon sheet. The equation 
which describes the gas density in the bomb as a function of time is: 

co p D q t 

2 Poo Y, ~2 ° 2 — exp ” 2 (8.1.3) 

n=1 P o " P 0 P °° + p oo^n £ 
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where : 


3 

p = density of gas in the bomb at any time, mg /cm 

3 

P q = unknown zero time density, mg/cm 

3 

P ot = equilibrium density in the bomb, mg/cm 
£ = half thickness of the sheet, cm 

2 

D = unknown diffusion coefficient, cm /sec 
t = time, sec 

= all non zero positive roots of the equation 

V- 

tan q = 

n p p 
o - °° 


(8.1.4) 


A generalized least squares method was applied to the pressure-time data 

for each test and the best estimate for p and D were obtained. 

o 

Figures 8.2 and 8.3 show data for two nitrogen tetroxide/TFE runs 
at different temperatures. The solid curve in each Figure was calculated 
from the best p^ and D while the circled points are the experimental 
data. 


8 . 2 Permeation Tests 

Two types of permeation tests were conducted during this study. 

The unidirectional tests were conducted by admitting the same gas or 
vapor to each side of a membrane but at approximately a 10 psi pressure 
difference. The counter diffusion tests were conducted with nitrogen 
or helium on one side and nitrogen tetroxide on the opposite side of 
the membrane at a near zero pressure differential. The method of 
operation and the method of data reduction differed for the two types 
of tests and thus are described separately. 
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Time, Seconds 

Figure 8.2 Adsorption of Nitrogen Tetroxide at 81°F 
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Time, Seconds 

Figure 8.3 Adsorption of Nitrogen Tetroxide at 163°F 
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8.2.1 Unidirectional Permeation 

The unidirectional permeation of a single gas, vapor or liquid 
through a membrane is widely measured and reported in the literature. 
Many designs of cells are described which produce scientific data, 
comparative data or go-no go results. The tests conducted as a part 
of this study were designed to provide scientific data and in particular 
the permeation constant defined by 


_ p A AP 

l t 


( 8 . 2 . 1 ) 


where : 


Q = quantity passed through the membrane in time (t) , CC(STP) or g 


P= the permeation constant, 


CC(STP) cm 


2 ° r 2 
sec cm atm sec cm atm 


_ J g_cm_ 


AP = pressure difference across the membrane, atm 

i = membrane thickness , cm 

2 

A = active area, cm 
t = time, sec 

A diagram of the unidirectional permeation apparatus is given in 
Figure 8.4. The fill lines and vacuum lines and associated valving were 
omitted for clarity. The test procedure consisted of the following 
steps : 

1. The bypass valve was opened and the cell and transducers were 
thoroughly degassed under vacuum. 

2. The test gas or vapor was introduced to the cell until the 
desired total pressure of the test was reached. In most tests 
one side of the cell was initially at vacuum and this step was 
eliminated. 

3. The bypass valve was closed and additional gas was introduced 
into the high pressure side (left side in the Figure) until 
the desired initial AP was established. 
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Figure 8.4 Unidirectional Permeation Apparatus 


8-9 




r tt 


4. The supply gas valve was closed and both the absolute pressure 
and the differential pressure were recorded as a function of 
time. Normally the tests were terminated when the differential 
pressure was one-half to three-fourths of its initial value. 

As in the case of the adsorption bomb tests described previously, 
the data reduction method for nitrogen tetroxide was different than the 
method for the permanent gases. When the change in differential pressure 
is small, it is possible to select an average AP and to obtain the 
permeation constant by an algebraic solution. However, when AP changes 
are large during the test and when a statistical treatment of the data 
is desired, the differential equation must be solved. For the permanent 
gases which obey the ideal gas low a closed form solution was obtained 


when: 

AP 

,o 


AP = 
R = 
M - 
T = 
A = 
£ = 
t - 


„ A ^ „ , „o RTAt /l , 1 x „ 

£n AP = In i P - -jjj- (y“ + y“) P 

H L 


differential pressure at any time, atm 

differential pressure at the start of the test, atm 

3 _ 

rti-/ cm atm 

gas constant, 82.054 - mo y— of 


( 8 . 2 . 2 ) 


gas molecular weight 
temperature, °K 

2 

active cross-sectional area, cm 
membrane thickness, cm 
time, sec 


V = volume on high pressure side, cm 
rl 

3 

V = volume on low pressure side, cm 

L 

P = permeability constant, ^ — 

sec cm atm 


The best value of p is determined by a least square treatment of the AP 
vs t data derived during the test. 
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For nitrogen tetroxide in which dissociation is a function of 
temperature and pressure, a closed form solution was not found. The 
differential equation which defines the permeation is the following: 


dm 

dt 


where : 


C 

C 


1 

2 


C 


3 


K 

P 

T 


H 


R 

o 

“l 

O 



(8.2.3) 


- K /8 
P 

4RT 

46.008 V u K 
H p 

= 4RT 

46.008 V T K 
L p 

= the equilibrium constant of dissociation 
= the absolute temperature, °K 


p N0 2 2 

Pn 2°4 


(8.2.4) 

(8.2.5) 

( 8 . 2 . 6 ) 


= the volume of the low pressure side of the permeation 
cell cm^ 

= the volume of the high pressure side of the permeation 

cell cm^ ^ 

r~ § cm atm 

= universal gas constant, 82.054 — 

g mole K 

= the zero time mass on the low pressure side of the 
membrane, g 

= the zero time mass on the high pressure side of the 
membrane, g 


Various values of P are assumed and the equation is numerically 
integrated. The value of P is selected which gives a least square 
error of pressure difference as a function of time calculated from 
the integrated values compared with the measured pressure differences. 
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8.2.2 Counter Diffusion 

The simultaneous counter diffusion of two gases through a membrane 
with a partial pressure difference but no total pressure difference is 
useful in determining whether interactions occur which affect permeation 
rates. The experimental apparatus which was used for this testing is 
shown in Figure 8.5. Again some accessory lines and valves were omitted 
to provide clarity. The test involved a continuous flow of the two 
gases from the supply tanks, through the cell and into or through the 
sample bottles. By careful adjustments of the several valves, the flow 
rate, the cell pressure and differential pressure were adjusted. After 
sufficient time was allowed for equilibrium to be established, a sample 
of each gas was taken for analysis. 

In all cases the partial pressure difference was equivalent to the 
absolute pressure of the cell and data reduction required only the 
algebraic solution of the standard permeation equation 


Q 


PA AP 
£ 


(3.2.7) 


where the terms were defined previously. 
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9.0 RESULTS 
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There are three basic classes of engineering data which were obtained 
during this study. Th?se classes are: primary (calculated directly from 

the experimental measurement) , secondary (calculated from the experi- 
mental data but requiring assumptions relative to mechanism) and derived 
(calculated from primary or secondary data) . It is evident that 
solubility is primary data, since the quantity of a gas which dissolves 
in a unit volume of a material is directly measurable. However, the 
Henry's Law constant (S) which is defined by 

q = SP (9.0.1) 


where : 


= the quantity dissolved per unit volume (usually 


P 

S 


CC(STP) _g_ 

3 3 

cm cm 

pressure, atm 
proportionality constant 


is secondary data which depends on the manner in which the solubilities 
follow the assumed law over the pressure range of interest. Similarly 
the diffusion coefficient computed from the transient adsorption data 
and the permeation constant obtained from the permeation cell tests 
represent secondary data. Examples of derived data include S, D or P 
calculated from the other two using the equation 


p = DS 


(9.0.2) 


where : - 

CC(STP) cm 

p= permeability constant, 2 

sec cm atm 


D = diffusion coefficient, cm /sec 

CC(STP) 

S = Henry Law constant g L 

atm cm of Teflon 





Thus, it should be noted that the engineering data required for the 
analysis of bladder leakage which was obtained as a result of the 
experimental effort is all of the secondary and derived types. In 
some cases the theory was found to be only marginally applicable and 
extrapolation of data outside of the range of measurement must be 
considered a dangerous practice. 

9.1 Solubility Data 

Henry's Law constants (S) were obtained from the bomb tests for 
nitrogen, helium and nitrogen tetroxide in polytetraf luoroethylene (TFE) 
and tetrafluoroethylene-hexaf luoropropylene copolymer (FEP) , Numerous 
attempts were made to obtain solubility data for hydrazine vapor but 
the gas phase decomposition was sufficiently rapid in the Cres steel 
(Types 316) apparatus to prevent meaningful data from being obtained. 

For both nitrogen (N^) and helium (He) , S was found to be inde- 
pendent of pressure over the range of pressure investigated and could 
be fit to the usual Arrhenius type equation for temperature with an 
apparent accuracy of about + 10%, Nitrogen tetroxide (NTO) data showed 
that S increased with increasing pressure but over the pressure range 
studied (i.e., about one-tenth to one-half the vapor pressure) an 
average value was selected which represented the data within about +25%. 

Figure 9.1 shows an example of the nitrogen tetroxide solubility 
data. The Henry's law constants which best represent these data give 
the straight lines which cross the data. Clearly the curved lines 
represent the data more accurately. 

The equations for these curves are: 


NTO /TFE 

, ,ry , n -3 1.36 

q = 3.42 x 10 p 

, 3 

g/cm 

(9.1.1) 

NTO /FEP 

-3 1 59 

q = 1.84 x 10 J p A,J * 

g/cm 3 

(9.1.2) 


These equations represent the experimental data within +2%. If S is 
defined as dq/dP, then 


NTO/TFE 

S = 4.65 x 10 3 P^* 3 ^ g/atm cm 3 

(9.1.3) 

NTO/ FEP 

_q o 59 1 

S = 2.93 x 10 P g/atm cm 

(9.1.4) 
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However, average values for S in the experimental range were used. These 
data were fit to an Arrhenius type equation to provide a good representa- 
tion of the effect of temperature. The experimental data are summarized 
in Table 9 in the form of Arrhenius type plots which represent the 
experimental data to approximately the indicated accuracy. 

Table 9-1 

Henry's Law Solubility Constants 

n 2 /tfe (+ 10%) 

S = 6.6 x 10~ 4 exp CC(STP) / atm cm 3 4-30 atm 

25-100°C 

n 2 /fep (+10%) 

_-i QOO 3 

S = 4.2 x 10 J exp CC(STP) / atm cm 8-30 atm 

25-100°C 

He/TFE (+10%) 

S = 2.5 x 10 -4 exp CC(STP)/atm cm 3 6-30 atm 

1 K. 

25-100°C 

He/FEP (+10%) 

S = 1.1 x 10 -3 exp CC(STP) /atm cm 3 8-25 atm 

25-100°C 

NTO/TFE (+25%) 

-7 3300 3 

S = 3.4 x 10 exp ^,o K ~ g/atm cm 0.1-0. 5 v.p. 

25-100°C 

NT0/FEP (+25%) 

S = 1.2 x 10 ^ exp 3 *°° g/atm cm 3 0.1-0. 5 v.p. 

25-100°C 
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9.2 Diffusion Coefficients 


m 


m 


Diffusion coefficients for gases and vapors in teflon were obtained 
both as secondary data and derived data. In this section only the 
secondary data, namely those obtained from the bomb tests, are discussed. 
Data were obtained for both nitrogen and nitrogen tetroxide with TFE 
and FEP. With helium the adsorption was so rapid, useful data points 
were not recorded and hydrazine vapor decomposed as discussed previously. 

Data for nitrogen diffusion were derived point by point and 
Figure 9.2 shows the adsorption curve for one test and the calculated 
point by point diffusion coefficients for this test and also for a 
duplicate test. The data were generally within +20% of the mean value 
for each temperature. 

Data for nitrogen tetroxide diffusion required a computer solution 
due to the dissociation reaction. These data are given in Tables 9-2 
and 9-3 for TFE and FEP respectively. There appears to be some pressure 
dependence for D, particularly in the NTO/TFE system, but the scatter 
is sufficient to limit this to an observation only. The high temperature 
tests generally contained an insufficient number of data points in each 
test to provide high confidence in any specific value of D. However, 
the 3o deviations are given to indicate the conformity of the data from 
each run to the best value of the calculated diffusion coefficient. 

The NTO/TFE data (except for the one point at 100°C) were fit to an 
Arrhenius equation which was then evaluated at the four test temperatures. 
The results given in Table 9-2 show that most data points are within 
+50% and the worst is about a factor of 3 from that best fit to the 
Arrhenius equation. The NTO/FEP data given in Tables 9-3 show less 
scatter. The lower scatter in the data is attributed to the lower value 
for D which permitted a larger number of data points to be taken during 
each test. For NTO/FEP the data generally is within +20% and the worst 
point is less than a factor of 2 from the best fit to the Arrhenius 
equation. 

The diffusion coefficient data obtained from the adsorption bomb 
test and the indicated accuracy are summarized in Table 9-4. 
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Elapsed Time, Minutes 
Figure 9.2 Adsorption of Nitrogen in 


Table 9-2 


Experimental D for NTO/TFE 


mp . °C 

P ave, atm 

D x 10^, 

cm^/sec 

Remarks 


27 

.15 

2.86 

+ 

.20 

15 

data 

points 

27 

.25 

4.77 

+ 

.63 

14 

data 

points 

28 

.3 

5.98 

+ 

.63 

12 

data 

points 

73 

1.4 

30.4 

+ 

4.0 

8 

data 

points 

73 

2.2 

33.7 

+ 

6.4 

5 

data 

points 

73 

2.9 

38.8 

+ 

1.4 

5 

data 

points 

74 

3.5 

66.3 

+ 

4.4 

5 

data 

points 

96 

5.8 

156 

+ 

6 

4 

data 

points 

96 

7.5 

113 

+ 

11 

4 

data 

points 

96 

8.7 

137 

+ 

14 

5 

data 

points 

96 

9.7 

186 

+ 

3 

4 

data 

points 

96 

8.0 

87 



2 

data 

points 

100 

3.3 

77 

+ 

13 

4 

data 

points 

100 

4.3 

63 

+ 

5 

4 

data 

points 

100 

5.1 

152 

+ 

560 

4 

data points 
(deleted) 

100 

5.8 

34 

+ 

2 

4 

data points 


Calculated 

from 

Arrhenius 

27 

4.5 x 10 8 

2, 

cm /sec 

73 

38 x 

00 

i 

o 

r*H 

2, 

cm /sec 

96 

90 x 

10- 8 

2, 

cm /sec 

100 

100 x 

oo 

1 

O 

rH 

2. 

cm /sec 
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Table 9-3 


Experimental D for NTO/FEP 



Temp . °C 

P ave , atm 

D x 10 

8 

5 

2. 

cm /sec 

Remarks 

— 

26 

0.3 

0.90 + 

.14 

Uncertainty in final 


28 

0.3 

0.81 

+ 

.13 

Temperature measure- 
ment, 9 data points 


29 

0.2 

0.97 

+ 

.07 

16 data points 


76 

2.3 

19.8 

+ 

2.4 

5 data points 


77 

4.4 

19.9 

+ 

3.3 

6 data points 

— 

77 

4.9 

36.1 

+ 

1.1 

5 data points 


78 

3.2 

17.9 

+ 

1.9 

7 data points 


78 

3.9 

21.2 

+ 

2.5 

5 data points 


92 

2.3 

40.3 

+ 

6.3 

6 data points 


92 

4.1 

41.6 

+ 

4.2 

5 data points 


93 

5.8 

53. 

+ 

15 

4 data points 


94 

7.2 

70. 

+ 

10 

3 data points 


95 

8.4 

54 

+ 

11 

6 data points 


Calculated from Arrhenius Equation 


28 

0.8 

X 

10 cm /sec 

77 

20 

X 

in -8 2, 

10 cm /sec 

93 

50 

X 

in” 8 2 / 

10 cm /sec 
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Table 9-4 


Diffusion Coefficients 


N 2 /TFE (+20%) 

D = 0.83 exp - cm^/sec 3-30 atm 

IK 

25 - 100 °C 


N 2 /FEP (+20%) 

-2 

D = 5 . 5 x 10 exp 

NTO/TFE (+50%) 

D = 0.40 exp - ^-P - 9 . cm^/sec 0.1 - 0.4 v.p 

25 - 100 °C 

NTO/FEP (+20%) 

D ® 56 exp - — ~ — cm /sec 0,2 - 0.4 v.p 

25 - 100 °C 


4000 2. 

^5^“ cm / sec 


10 - 25 atm 
25 - 100 °C 
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9. 3 Permeability Constants 


Unidirectional permeation measurements were made with helium, 
nitrogen and nitrogen tetroxide through TFE and FEP membranes. The 
results of these measurements are given in Table 9-5. The helium data 
gave permeability constants that were independent of total pressure 
and differential pressure and were a function of temperature only. 

Except for two tests, nitrogen gave similar results. 

The nitrogen tetroxide data presents a more complicated picture. 

The calculated permeability constants in both TFE and FEP show no 
clear dependence on temperature but are strongly influenced by the total 
pressure. The pressure dependence was computed by least squares and 
gave the following relationships: 

NTO/TFE px 10 10 = 5 + 8 P atm j a = 23 x 10 -10 

sec atm cm 

NTO/TFE px 10 10 = 3.7 + 2.3 P atm j a = 1.8 x 10” 10 

sec atm cm 

Clearly in the case of NTO/TFE, the data have excessive scatter as 
indicated by the high value of a and the equation cannot be used with any 
confidence. Additional testing at high pressure is required to define 
more accurately the pressure dependence. 

Counter flow permeation measurements were made with helium and 
nitrogen tetroxide through TFE and FEP. These data are presented in 
Table 9-6 in the same format and order as the unidirectional data of 
Table 9-5, The permeability data for all systems including nitrogen 
tetroxide appear to have a tendency to increase with increasing tempera- 
ture, There are a sufficient number of exceptions to limit this to an 
observation only. 

The counter diffusion data generally shows that helium permeates 
more rapidly by a factor of 2 or 3, nitrogen tetroxide permeates more 
slowly at low pressure by a factor of 2 to 5 and nitrogen varies in 
both directions within a factor of about 2. 
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Table 9-5 


Experimental Unidirectional Permeation 


Temp , °C 

P ave, atm 

AP , atm 

CC(STP) cm 
P, 2 

sec cm atm 

He/TFE 




23 

0.6 

0.6 

3.4 + 0.01 x 10 

23 

3.3 

0.6 

3.4 + 0.1 

65 

0.6 

0.6 

8.1 + 0.4 

74 

18.3 

0.6 

7.9 + 0.1 

86 

1.4 

0.8 

8.2 + 0.6 

He/FEP 




21 

17.5 

0.2 

2.5 + 0.2 x 10' 

39 

0.3 

0.2 

2.8 + 0.02 

39 

0.7 

0.7 

2.7 + 0.2 

60 

18.9 

0.7 

4.3 + 0.1 

65 

0.6 

0.6 

4.8 +0.5 

90 

0.8 

0.8 

9.8 + 0.5 

n 2 /tfe 




29 

4.7 

0.7 

1.5 + 0.2 x 10" 

38 

0.6 

0.6 

1.7 + 0.7 

69 

0.6 

0.6 

6.1 + 0.1 

69 

5.8 

0.7 

9.8 + 0.7 

77 

19.3 

0.7 

8.5 + 0.5 

90 

0.7 

0.7 

1.4 + 0.1 


9-11 


Table 9-5 (cont.) 


Temp , °C 

P ave, atm 

AP, atm 

CC(STP)cm 

P » 2 

sec cm atm 

N 2 /FEP 

35 

5.6 

0.7 

2.2 + 0.5 x 10' 

37 

0.3 

0.2 

0.8 + 0.3 

38 

17.6 

0.6 

11.2 + 0.3 

39 

5.1 

0.7 

1.7 + 0.8 

40 

1.1 

0.7 

0.5 + 0.1 

62 

0.7 

0.7 

3.2 + 0.3 

62 

19.3 

0.7 

4.7 + 0.4 

85 

17.8 

0.7 

7.1 + 0.8 

91 

0.5 

0.5 

10.3 + 0.2 

R cm 

p> 2 

sec cm atm 

NTO/TFE 

26 

0.6 

0.6 

5.1 + 0.5 x 10' 

60 

0.7 

0.6 

7.9 + 0.5 

61 

4.9 

0.5 

84. + 1 

64 

0.5 

0.5 

4.7 + 0.1 

85 

0.6 

0.6 

6.2 + 0.4 

87 

7.9 

0.7 

45. + 8 

NTO/FEP 

25 

0.5 

0.5 

5.1 + 0.8 x 10' 

64 

5.2 

0.6 

18 + 1.5 

65 

0.5 

0.5 

3.4 + 0.6 

84 

0.5 

0.5 

4.9 + 0.1 

84 

8.1 

0.5 

21 + 0.3 
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Table 9-6 


Experimental Counter Flow Permeation 


Temp , °C 

P ave, atm 

AP, atm 

CC(STP) cm 
P> 2 

He/TFE 

25 

0.3 

0.6 

sec cm atm 
9.2 x 10~ 7 

62 

1.0 

2.0 

23. 

87 

0.6 

1.3 

26. 

89 

1.0 

2.0 

27. 

He/FEP 

25 

0.3 

0.7 

8.7 x 10“ 7 

30 

0.3 

0.7 

1.3 

88 

0.3 

0.7 

8.9 

88 

0.6 

1.3 

7.6 

88 

1.0 

2.1 

16. 

n 2 /tfe 

25 

0.3 

0.7 

12 x 10 -8 

87 

0.6 

1.3 

(questionable 

analys: 

2.3 

87 

1.0 

2.0 

3.8 

n 2 /fep 

88 

1.0 

2.0 

13. x 10“ 8 
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Table 9-6 (cont.) 


-T -' 

Temp , °C 

P ave, atm 

AP , atm 

P> g cm 





2 

sec cm atm 


NTO/TFE 




? "* 




-10 

r 

25 

0.4 

0.8 

0,29 x 10 


25 

0.3 

0.7 

0.42 


62 

1.0 

2.1 

2.4 


87 

0.6 

1.3 

0.96 


87 

1.0 

2.0 

1.25 

— 

87 

0.6 

1.3 

1.67 


89 

1.0 

2.0 

3.1 


NTO/FEP 




— 




„ , „-10 


25 

0.3 

0.7 

0.18 x 10 


30 

0.3 

0.7 

< 0.1 

^ . 

88 

0.3 

0.7 

0.78 


88 

0.7 

1.3 

< 0.09 

fc : 

88 

1.0 

2.1 

1.03 


88 

1.0 

2.0 

1.35 
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A review of all the experimental and analytical steps involved in 
making the permeability measurements showed no clear reason for the 
variability of the data. The least squares determination of the uni- 
directional permeation data for nitrogen tetroxide reported in Table 9-5 
(see Section 8.2) showed that p is not constant but decreases during 
each run. The reported values are the best constant which fits the 
data over the pressure range for each test. It was estimated from the 
apparent pressure sensitivity that p would not vary by more than about 
20% during any test. The effect on the tabulated data would probably 
be small since an average pressure, rather than the final, pressure 
was reported. Table 9-7 summarizes the permeability results and 
indicates the range of conditions over which the data were obtained. 
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Table 9-7 


Permeation Constants 


He/TFE 

n nn in -A 1650 CC(STP)cm 

p = 0.88 x 10 exp - ^ok — 2 

sec cm atm 


_ -7 CC(STP)cm 

p = 9 to 27 x 10 1 — z 

c 2 

sec cm atm 


He/FEP 

. „-4 1950 CC(STP) cm 

p = 1.5 x 10 exp - ^ — 2 

sec cm atm 


. .. . . , . «-7 CC(STP) cm 

= 1.3 to 16 x 10 * — 2 

sec cm atm 


N 2 /TFE 


„ _ -2 4300 CC(STP)cm 

p = 2.0 x 10 exp - — 2 

sec cm atm 


= about 3 x 10 


-8 CC(STP) cm 
2 . 

sec cm atm 


n 2 /fep 


, „-2 4500 CC(STP)cm 

P u = 2.4 x 10 exp - — 1 — 2 

sec cm atm 


. . ... ,„-8 CC(STP)cm 

p - about 13 x 10 7) 

c 2 

sec cm atm 


NT0/TFE 


Pu 


about 6 x 10 


-10 g cm 


= 0.3 to 3 x 10 


sec cm atm 

-10 g cm 

2 _ 

sec cm atm 


0.5 - 18 atm 
20 - 90° C 

0.3-1 atm 
25 - 90°C 


0.3 - 19 atm 
20 - 90°C 

0.3-1 atm 
25 - 90°C 

0.6 - 19 atm 
30 - 90°C 

0.6-1 atm 
85°C 

0.3 - 19 atm 
35 - 90°C 

1 atm 
90°C 

0.6 atm 
25 - 90°C 

0.3-1 atm 
25 - 90°C 
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Table 9-7 (Cont.) 


NTO/FEP 

i_ ^ ^ *1 ^“10 r cm 

n = about 5 X 10 “ — — 

nj 2 

sec cm atm 


0.5 atm 
25 - 90°C 


= 0.2 to 1.4 x 10 


-10 


cm 


sec cm atm 


0.3-1 atm 
25 - 90°C 


p u = permeability constant 
p c « permeability constant 


from unidirectional tests 
from counter flow tests 
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9.4 Summary of Results 


h’H 






The preceding sections presented the experimental results separately 
for each type of test conducted during the experimental program. In this 
section the complete data generated by the experimental program is inter- 
compared and then compared with the available literature. Since only 
permeability constants were found in the literature, it is most convenient 
to examine all data in terms of the permeability constant. This constant 
for many flow geometries is directly proportional to the leakage rate. 


9.4.1 Helium Permeation 


Helium permeation data are available from the experimental program 
for both counter diffusion and unidirectional permeation through TFE and 
FEP. Since diffusion coefficients for helium could not be measured during 
the adsorption tests, there is no permeation data derived from the SD 
product . 


Stern et al (Reference 18) give an Arrhenius type equation for 
helium permeation in FEP which applies over the measured range of 0.1 to 
4 atm and 25 to 100°C. Their equation reduces to the following: 


"He 


_3 

1.4 x 10 exp 


2400 CC(STP)cm 

T 2 

sec atm cm 


In an earlier article, Stem et al (Reference 19) gave data for helium 
permeation in TFE in graphical form over the same temperature range. 
Graphical data from duPont was also reported by Redel (Reference 22) in 
the same temperature range. 

The only other data found for helium permeation of Teflon is from 
Rocketdyne (Reference 21) where data is given for a mixed sheet of TFE 
and FEP and for TFE sheet containing aluminum flakes. This cannot be 
readily compared with data for pure materials. 

All the available permeation data are summarized in Figure 9.3. 

The literature values are represented by the solid lines and the least 
square fits of the unidirectional permeation data are represented by the 
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r -i 

|^3 


i fci 


j 


dashed lines. The actual unidirectional data are shown as open circles 
and triangles and the counter flow data by the corresponding filled 
symbols. The literature data is in good agreement with the data obtained 
during this study. 


9.4.2 Nitrogen Permeation 

Nitrogen permeation data are available from the experimental 
program for both counter diffusion and unidirectional permeation through 
TFE and FEP. It also was possible to compute the S*D product for both 
TFE and FEP from the data measured in the adsorption bomb tests, namely: 


n 2 /tfe 

n 2 /fep 


S'D = 5.5 x 10 ^ exp 
-4 

S*D = 2.3 x 10 exp 


3000 CC(STP) cm 

T . 2 

sec atm cm 

3100 CC(STP) cm 

T 2 

sec atm cm 


Stem et al (Reference 18) give an Arrhenius type equation for nitrogen 
permeation in FEP which is applicable at pressures from 0.1 to 4 atm at 
temperatures from 60° to 100°C. Below 60°C it is applicable at low 
pressures and higher, at least as high as 4 atmospheres if sufficient 
time is allowed for equilibrium to be attained. At room temperature 
the times required for equilibrium appears to be of the order of a few 
hours. Their equation reduces to the following: 

_ _ ,„-3 3600 CC(STP)cm 

p = 3.2 x 10 exp - — * L j 

sec atm cm 

In an earlier article Stem et al (Reference 19) gave data for nitrogen 
permeation in TFE in graphical form over the range 25 to 90°C. No 
mention was made of pressure dependence at low temperature with TFE. 

—8 

Lebovits (Reference 17) gave a value of 1.6 x 10 
2 

CC(STP) cm/sec atm cm for the permeation constant of nitrogen through FEP 

at 20 to 30°C and Rocketdyne reported a value of 4.6 x 10 ^ 

2 

CC(STP)cm/sec atm cm at 21 to 2A°C. Graphical data from duPont for 
nitrogen permeation of FEP over the temperature range of 0 to 100°C was 
reported by Redel (Reference 22) . 
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B 




m 


! 
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These data are summarized in Figure 9-4. The literature data are 
represented by solid lines, the least square fits of the unidirectional 
permeation data obtaired during this study are represented by dashed lines 
and the SD products are shown by dotted lines. The unidirectional and 
counter flow data points are shown in the figure as open and filled symbols 
respectively. The literature data and those data obtained during this 
study are in good agreement. Tn view of the discussion of Stern (Refer- 
ence 18) of the long times required for ^/FEP to reach equilibrium at 
pressure as high as 4 atmospheres and at temperatures below 60 C, the 
unidirectional data at 35 to 40°C should be noted. These data range from 
5 x 10 to 1.1 x 10 CC(STP) cm/sec atm cm for average pressures ranging 
from 0.3 to 18 atmospheres. The highest value taken at 18 atm was not 
included in the curve fit which is believed to represent the equilibrium 
values. Considering the probable precision of the experiments, it is 
likely that this point represents a nonequilibrium value of the perme- 
ability constant. 


9.4.3 Nitrogen Tetroxide Permeation 


Nitrogen tetroxide permeation data were obtained in the experimental 
program for both unidirectional and counter flow through TFE and FEP. 

Since the adsorption tests provided both solubilities and diffusion 
coefficients, it was also possible to obtain the corresponding SD product 
namely : 


NTO/TFE 

NTO/TFE 


S*D = 1.4 x 10 ^ e xp 
S • D = 6.7 x 10 ^ exp 


1500 g cm 

T . 2 

sec atm cm 

3200 g cm 

T 2 

sec atm cm 


The literature also provided some permeation data but these data 
were available at only one or two temperatures and in a form which requires 
assuming the differential pressure to derive the permeability constant. 

JPL (Reference 20) reported the results of a test at 20°C through a 10 mil 
teflon film (assumed to be TFE) which showed an equilibrium leakage rate 
of 5.6 X 10" g/ hr in . From the nature of the experiment, it was 
possible to estimate an equilibrium permeability constant of 
6.5 x 10 g cm/sec atm cm . 
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Figure 9.4 Experimental Permeation of Nitrogen 
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Red el (Reference 22) reported leakage rates of 1.7 to 4.7 x 10 

g/hr in^ through single 2 mil sheets of FEP at 24 C and 7 x 10 g/hr in 

through double 2 mil sheets of FEP. These leakage rates permitted 

_9 

permeability constants to be calculated ranging from 3.0 x 10 to 

—8 2 
1.0 x 10 g cm/sec atm cm . 

Rocketdyne reported the results of a large number of tests using 
several test procedures and several pure and composite samples. The data 
for pure materials are given in Table 9—8 together with the calculated 
permeability constants. 


Table 9-8 

NTO Permeability (Reference 21) 


Material 

u 

o 

H 

Q/A, 

lbs /hr in^ 

p calc. 

JL 

cm/sec ; 

10 

mil 

FEP 

74 

6.76 

x 10 -5 

4. 

.0 

X 

10" 

-9 

10 

mil 

FEP 

74 

6.13 

x 10 -5 

3. 

,7 

X 

10" 

■9 

10 

mil 

FEP 

74 

13.7 x 10~ 5 

8. 

.2 

X 

10" 

-9 

10 

mil 

FEP 

74 

13.6 x 10 -5 

8. 

,2 

X 

10" 

-9 

15 

mil 

TFE 

74 

4.67 

x 10" 5 

4, 

.2 

X 

10" 

-9 

15 

mil 

TFE 

25 

0.29 

x 10~ 5 

2. 

,2 

X 

10" 

■9 

10 

mil 

FEP 

21 

0.43 

x 10" 5 

8, 

.0 

X 

10" 

-10 


All nitrogen permeability data obtained to the experimental program 
and calculated from the literature are given in Figure 9.5. The uni- 
directional and counter flow data were discussed in detail in the previous 
section but it is more evident from the figure that unidirectional tests 
gave permeation constants up to an order of magnitude higher than the 
counter flow tests. The S*D product, shown by the dotted lines is in 
better agreement with the unidirectional data as would be expected from 
the nature of the adsorption test. Most of the literature values gave 
results higher than all but the high pressure unidirectional tests of 
this study. It was not always possible to determine whether the litera- 
ture test represented counter flow or unidirectional. The JPL data and 
apparently the Redel data and some of the Rocketdyne data are counter flow 
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X TFE^-.TFE (S*D) 


▲ 




V FEP (s-d) 


A 

▲ 


Nitrogen Tetroxide Permeation 

0 FEP unidirectional 
• FEP counter flow 

A TFE unidirectional 
▲ TFE counter flow 
X TFE finat rate Ref. 20 

1 FEP calculated from Ref. 21 
X TFE calculated from Ref. 21 
^ FEP calculated from Ref. 22 





V 


fej 


results. The high permeation rate data of Rocketdyne at 74°C are 
apparently counter flow but the remainder may be unidirectional. In all 
cases, the literature data seems to have liquid nitrogen tetroxide rather 
than vapor in contact with the membrane and therefore the test pressures 
at each temperature were above those of this study. 

9.5 Leakage Through Holes 

A final part of the test program involved counter flow tests with 
a laminated sample containing a selected hole area. The sample was made 
from 1 mil aLuminum foil containing 100 ten mil diameter holes spaced on 
0.1 inch centers of a 10 by 10 square pattern. To this foil 10 mils of 
TFE and FEP were applied to opposite faces. The sample was placed in the 
permeation cell and run in the counter diffusion mode with nitrogen 
tetroxide on the TFE side of the sample and helium on the FEP side. 

During the first of three tests, the sample was evacuated at 187°F 
and when thoroughly degassed nitrogen tetroxide vapor and helium gas 
were introduced to a total pressure of 145 psia on each side of the 
sample. The cell was closed off and after 2 hours the gases on each side 
of the membrane were withdrawn for analysis. This test gave a helium 
leakage rate of 1.5 x 10 CC(STP)/sec, but an undetectably low concentra- 
tion of NO^. Equilibrium clearly was not attained during this time 
period. 

The two later tests were conducted with a modified procedure. After 
evacuating the sample, the nitrogen tetroxide and helium were slowly 
flowed past the membrane at approximately the pressure to be used during 
the test. After a period of several hours, the cell was closed off as 
in the first test and allowed to remain for the specified test duration. 
The gas from each side was withdrawn for analysis as in the first test. 
Helium analyses in nitrogen tetroxide were by a gas chromatographic 
method with nitrogen as a carrier gas and nitrogen tetroxide analyses in 
helium were by an infrared technique. 
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The results of these tests are as follows: 


Test 2 Test 3 


Temperature, °C 

86 

93 


Pressure, atm 

9.9 

4.9 


Flow time, hours 

6 

2 


Test time, hours 

Helium leak rate, CC(STP)/sec 

64.5 

8.6 x 10" 8 

44 

1.6 x 

IQ" 5 

Nitrogen tetroxide leak 
rate g/sec 

3.1 x 10 -8 

1.6 x 

io -8 


The membrane was not inspected between tests, but at the conclusion 
of the tests the laminated sample was removed and examined. Several 
changes had occurred. The membrane was supported on the TFE side by a 
screen and on the FEP side by a porous metal sheet so that gross move- 
ment due to momentary pressure differentials was restrained. Neverthe- 
less, at some time during the test series sufficient movement occurred 
to tear the aluminum foil near the flange. This may have occurred when 
flow rates were being adjusted at the start of Test 2 or Test 3, or 
possibly when the gas samples were taken at the end of any of the three 
tests. At any rate, the time at which this additional leakage path 
through the barrier occurred cannot be determined. 

Another change also was evident. There was complete separation of 

2 

the TFE and aluminum over the full 45.5 cm area of the sample. This 
separation could have occurred at any time during the test series, but 
it is likely to have occurred during the initial test either during 
evacuation or during depressurization as the gas samples were withdrawn. 
Thus, very likely Tests 2 and 3 were made with only the FEP bonded to 
the foil. The sample as received had a weak bond at the TFE-foil inter- 
face. As in all the previous tests involving FEP membranes, the FEP 
was considerably lighter in color after exposure to nitrogen tetroxide 
but no significant effect was attributed to the bleaching. Earlier 
laminated samples which had insufficient hole areas were tested in the 
adsorption bomb with similar results relative to bleaching and delamina- 
tion. It appears that improvements in bonding between TFE and aluminum 
foil are necessary. 
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Appendix A 


Consider a laminated bladder structure separating two infinite half 
spaces. The upper half space contains gas, the lower liquid, and the 
laminated structure consists of two layers of material permeable to the 
pressurant gas separated by an impermeable barrier which has one circu- 
lar hole in it. We wish to determine the leakage rate of the gas 
through the bladder and into the liquid taking into account the effect 
of the gas and liquid spaces. Let p^ be the partial pressure of the 
pressurant gas in the gas mixture, c^ the concentration (in mass fraction) 
of the gas in the first layer of Teflon, c 2 that in the second layer, 
and c^ its concentration in the liquid propellant. Then from figure 
(A. 1) i x 


0 , x > h 

0 , 0 < x < h 

0 ,-b < x < 0 

0 , x < -b 



t r 


Gas 

Barrier 

Material 1 


Material 2 


Liquid 


x = +h 
x - 0 


Figure A.l 


Instead of p^ we solve for = P Q - where p^ is the total 
system pressure. Since p^ goes to p at distances far removed from the 
opening, vanishes for large (x +r ) ^ • In cylindrical coordinates* 


V 2 *. + 0 

Y .2 r 3r „ 2 

3r 3x 

The Hankel Transform of zero order of eq. (A.l) gives 
? 

, /1 \ n 


L ^- - k^ »J>(k,x) = 0 


(A.l) 


(A. 2) 


as the equation satisfied by the Hankel Transform of Kr,x) 
solution which vanishes for large Vx 2 + r^ is 


*We assume a circular shaped hole for convenience. The results are 
however valid for holes of any shape. 


A-l 


(A. 3) 


i|)(k,x) = A^(k)e 


-kx 


where, of course, A^(k) is an unknown function of k. Forming the 
Hankel inversion and using the definition of i|»(r,x) gives 


p^r.x) = p Q - 


— W V 

kA, (k)e J (kr)dk , x > h 
1 o 


(A. 4) 


For the Hankel Transform of c^(r,x) we find in the same manner 


c^(k,x) = A 2 (k) cosh kx + A^k) sinh kx 


(A. 5) 


with inverse. 


c ^(r ,x) 


kA 2 (k) cosh kx J Q (kr) dk 


+ 


kA^(k) sinh kx J Q (kr) dk 
' o 


0 < x < h 


(A. 6) 


The boundary conditions which must be satisfied at x = h are 


and 


c 1 (r,h) - K 1 p 1 (r,h) 

3 Cl (r,h) M D ^ 

n D = = — - Q - — 

P 1 1 Bx RT 9x 


(A. 7) 
(A, 8) 


In the first equation is the solubility of pressurant gas in material 1 
while in the second is density of material 1 and is the diffusivity 
of pressurant gas through material 1. Also M is the molecular weight 
of the pressurant, the binary diffusion coefficient for the mixture 
of pressurant gas and propellant vapor, R is the gas constant and T 
the gas temperature. Applying eqs. (A. 7) and (A. 8) to equations (A. 4) 
and (A. 6) we obtain from (A. 7) 


A- 2 


fOO fOO 

I kA 2 (k) cosh kh J Q (kr) dk + I kA 3 (k) sinh kh J Q (kr) dk 

'o' 

f 00 

= K . p - k, kA.(k)e“ kh J (kh) dk (A. 9) 

1 o 111 o 

* n 


and from (A. 8) 


PD, 


k 2 A„(k) sinh kh J (kr) dk + k A (k) cosh kh J (kr) dk 


L'O 


M D 
_S_S. 
RT 


kV (k)e kh J (kr) dk 
1 o 


(A. 10) 


These equations hold for all r. We may, therefore, write from the 
second equation 


A^(k) sinh kh + A^(k) cosh kh = 


M D 
- & Sl. 


P^RT ‘ 1 


A, (k)e 


-kh 


(A. 11) 


Using the representation of the delta function 

r r J o (kr)dr = ^ (A. 12) 

'o 

we obtain, after multiplying equation (A. 9) by r J Q (k,r) and integrating 
over all r, 

A 2 (k) cosh kh + A^(k) sinh kh + <]A^(k)e = K-jP Q £ ^ (A. 13) 

From eqs . (A. 11) and (A. 13) we eliminate A^k) and solve for A^(k) 
in terms of A 2 (k) . Substituting the value of A^(k) in terms of A 2 (k) 
in eq. (A. 5) we obtain 

. . *lP 0 sinh kx 6 (k) 

c l^ ,x ' ~ tanh kh + 8 cosh kh k 

+ V k > [ COSh kx ' tanh^kh"+ F slnh kx ] (A ' U> 
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e 



EJf 



; *2 




where 8 = k.o,D,RT/M D . If 6 is very small, that is, zero, eq. (A. 14) 
11 1 g g 


reduces to 


sinh kx 6(k) sinh k( h-x) 

c l (k ’ x) = K l P o sinh "kh k + sinh kh A 2 (k) (A ' 15) 


which, upon inversion gives 


x 


c l (r * x) = K l P o h + J 


kA 0 (k) Slnl ? — J (kr) dk (A. 16) 

2 sinh kh o 


Putting 


A 2 (k) = K l P o 


- A 2 '(k) 


gives finally 


c 1 (r,x) = k 1 p o 


kA 2 ’(k) 
' o 


sinh k(h-x) 
sinh kh 


J (kr) dk 
o 


(A. 17) 


This is the same equation for A ^ as was found in section 2,1 for the 
case of constant concentration, c q = k^p q , on tbe boundary x = h. We 
conclude that when 3=0 the gas phase can be ignored and the problem 
reduces to that with a constant concentration on the gas phase boundary 
of the bladder. 


Let us examine the order of magnitude of 3 under typical conditions. 
For Nitrogen gas the measured value of the solubility is of the order 
of 10 ^ [(gms N 0 )/(cm^) ] /[(gms Teflon) /(cm'*) } atm \ At 300 K 

z n _ £ 

- RT/M = (82.05) (300)/28 = 0.88 x 10 while v 10 since the 

measured value of D 1 is about 10 ^ and for binary gas solutions D 
is always about 10 . The density of Teflon is of the order of 1.0 so 


6 = 


k i p i d i rt 

M~D 


is of the order of 10 ^ . Thus 6 is expected to be very small. 
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The leakage rate must, for a physically meaningful solution, be 
a continuous function of the parameters, for example 6, upon which 
it depends. Therefore the leakage rate for very small values of 8 
will differ only slightly from its value for 6=0. As we have shown 
above 6=0 corresponds to the assumption of a constant gas concentra- 
tion on the interface between the bladder and the pressurant 8 as space. 
Therefore, the leakage rate through the structure can be calculated to 
good accuracy by ignoring diffusion in the gas space. 

A similar analysis for the interface between bladder and liquid 

space gives the result that the concentration on this interface can 

be assumed to be zero provided (P iPn^n) >> Here is 

the density of the liquid propellant, D £ is the diffusivity of pressurant 

gas in liquid propellant, is the solubility of the gas in material 2, 

P_ is the density of material 2 and k is the solubility of the gas in 
2 -5 2 
liquid propellant. In this case P £ e P 2 , * 10 cm /sec (a best 

-4 

estimate from the literature, see reference 1) ^ 10 as above, 

D„ ^ 10 ^ as above, and k p ^ 10 ^ according to the data of Chang and 

2 / I rt\ 2 

Gokcen^ . Thus the magnitude of the ratio is about 10 . It appears 

that neglecting diffusion through the liquid is not so good an approxima- 

2 

tion as neglecting it through the gas. However, 10 is sufficiently 
large compared to 1.0 for us to expect that calculations based on 
assuming the concentration vanishes at the liquid-bladder interface 
will give a good engineering approximation to the leak rate. 
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Appendix B 


In this section the integrals necessary for the calculation of the 
leakage rate through rectangular slits (Section 2.2) will be evaluated. 
From eq. (2.2.18) the integral is 


r 

J n 


dk 


tanh kh J (ka) sin ka — r = F 
o t ^ 

o k 


(B.l) 


2 2 

Let ka = u; dk/k - a(du/u ) ; kh » (h/a)u = au where a = h/a. 
(B.l) is equivalent to 

r , v du 

F(a) = a tanh au J (u) sin u — j 
^ o 0 u 


Then 


(B . 2) 


To evaluate the integral it is necessary to represent tanh au by its 
expansion in partial fractions^^ . 


tanh au 
u 


i T I 

01 n=o u 2 + (n+i) 2 (ir 2 /a 2 ) 


(B. 3) 


Using (B.3) in (B.2) gives, after interchanging the order of summation 
and integration 

J (u) sin u 


r , , !»vr v u > sln u du 

™ 2 7 2 — — 
n=o J o (u + Y n ) 


(B . 4) 


1 7T 

where y = (n + x) - . Since the integrand is an even function of u, we 
? n 2 a 

can write eq. (B.4) as 


F(a) = — / . f (a) 

a » — ^ n 

n=o 


(B.5) 


where 


f n (a) 


/4<» J (u) sin u m 
' o du 

,2 2v u 

( u + Y n > 


B-l 


v-i 


s 

to 


jEI 
i r I 

w 


E 

B 


I i 



j © 


! pi 

S B 


To evaluate the f (a) consider the complex integral 
n 


J (z)e 
o 


iz 


dz 

c (z2 +Y n 2) Z 


where z = u 4* iv, and the contour C is a semi-circle in the upper half 
complex plane indented at the origin as shown. 



By the theory of residues we have 


f V z)elz dz = r e V u) elu du f J o (z)elz 

L z 2 + 2 Z J-R u(u 2 + y 2 ) '^(z 2 + y 2 ) 

n n 1 n 


dz 

z 


R J (u) e 
o 


(u 


iu t / \ iz 

e . /• J (z)e 

du + _2 

+ Y n 2 ) u ■'r^z 2 + ^ 2 ) 


dz 

Z 


= 2ui T p (B.6) 

n 

where p denotes the n^ residue of the integrand within the contour, 
n 

On T (the large semi-circle) from the asymptotic expansion of 

t / 7(11) U 

J (z) we have 

o 




B-2 


J (z) ~ z e ^ z cos(z - tt/4) 

o 


-1/2 iz 
z e 


f .«-? + 

l ; 


Since on z = R(cos 0 + i sin 0) # 


iz, , N 1 , -2 z sinfl . 

e J (z) < — TTo~ (e 1 1 +1) 

° - I | 


e “2 1 2 I sin 0 x 

u 1/2 i + “u I7y i 


Now sin 0 is positive in the interval 0 < 0 < tt so that 
e ^ as | z | becomes large everywhere on the contour except 

possibly at 0 = 0 and 0 = tt. Thus everywhere on F ^ 


iz T , x A 

|e yz)^-^ 

R 


1/2 1/2 

where A is a finite constant and R = z . Then 


elZj 0 (2> A 

z(z 2 + 2 > V /2 

n 


Since the length of the contour is ttR, the integral over V ^ vanishes 
as R — >a like 1/R^^. 
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I i 

u 


u 


The only singularity of the integrand in the upper half plane is at 
a simple pole at z = iy . The residue of the integrand there is 


J (iy )e Y n 
o n 


I (Y) e ~ Yn 
o n 


(B . 9) 


2y, 


n 


where I (y ) = J (iy ). Using the above results in eq. (B.6) gives the 
on on 

following preliminary results. 


f-e T , x iu 

J o (u)e du . 

V z)eiz 


rR T , x iu 
J o (u)e 

du 

2 2 u + 

J-E U +1r „ J 

r z 2 + y 2 
r l n 

— + 
T 

Z 

4 

2 2 
u + y 

n 

e 

u 


= -2iri 


1 (y )e~ Yn 

o n 

2y 2 


(B.10) 


_ ^ i0 

On r i put z = ee 

J (z)e lz , 
o dz 

2 , 2 z 

r z + y 
r l n 


j A 

dz * ee i d0 to get 


J (ee 10 ) inJ (0) 

° i dO = - — — 


iir 


2 10 2 
_ e e 4* y 
0 =tt n 


'n 


t 


for small e, since J^CO) = 1* Letting e 0 and R > 00 , we obtain the 
result 


i m 


Jo (u)e 1 " du Jji lj_ j j ) -Y 

2 2 u 2 2 o' V 

u + y Y Y„ 

oo 'n n n 


(B.ll) 


Where the P denotes that the principal value of the integral is 
intended. On expanding the exponential in the integrand eq. (B.ll) 
becomes 
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\ 

u 


+ a. 1 + ~ 

cos u J (u) , I J (u) sin u , 

o du _o du 

2 2 u J 2 2 u 

u + T ■'-*°u+Y 


*72-72 W"~ T *J 

Y n 


T°o 

( cos u 

"77 


cos u J (u) J 

2 ^ = 0 

2 . 2 u 


f n (a) = 


J (u) sin u , 
o du 

2 2 u 

u + y 

n 


~2 


where the P has been omitted from the last Integral since it exists in 
the ordinary sense. 


From eq. (B.5) 


£ -Vv £ 

n=o y n=o 

*n 


I(y_)e Y n 
o n 


The first term is summable. From 


tanh au _ _2 ^ 1 

u ~ a 2 2 

n=o u +y 
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we have 


(B.16) 


llm tal,h ° u - a - 1 

U— K) U 3 


E 

n=o y 


Hence 


. aira air V'* 

F(a) = -5- - — Z-» 


11=0 


V T n )e Yn 


(B.17) 


2 2 

2 . . U jC 


Substituting the value of o = h/a and = (n + -y) 8*- ves 


F(a) 


«° T r/ . h irai -(n^lT 

irh h V* I o [(P 2 } h ]e 2 h 


(B.18) 


n=o 


(n + y) 


Using this result in eq. (2.2.19) and multiplying numerator and 
denominator by 2/v we obtain the result expressed in eq. (2.2.20). 

The function denoted by G(a) in section 2.2 is, of course, the second 
term on the right of eq. (B.18) multiplied by 2/irh. 
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Appendix C 

UNIDIRECTIONAL PERMEATION OF NITROGEN TETROXIDE 


fj 




IS 



i m 


The permeability constant p is defined by the differential 
equation 

4 s = 4 p ap (c.D 

dt Z 

where 

m is the net mass transferred from time zero (a function of time) 

4^- is the net mass transfer rate (a function of time) 
dt 

A is the area of the membrane (constant) 

H is the membrane thickness (constant) 

AP is the pressure differential across the membrane (V - P L ) 

(a function of time) 


The pressure on each side of the membrane is a function of the 
volume, the temperature and the total number of mole of gas* The 
volume and temperature are fixed by the experimental conditions* The 
number of moles on each side however, varies with the mass. Using the 
high pressure side for example the relationship between pressure and 
mass can be developed. Defining: 

T temperature, °K 

3 

V„ volume on high pressure side, cm 

H 

R gas constant, 82.054 atm cra^/g mole °K 

P u pressure, atm 

n^ g moles NgO^ 

n£ g moles N0£ 

n^ total moles (n^ + nj) 

K equilibrium constant for reaction NjO^^i 2 N 02 * atm 
(defined by T) 

hl, total mass of gas on the high pressure side, g 


C-l 


Perfect gas law 


(C. 2) 


P H V H n T RT 


Total Mass 


"h 


92 ,016n, + 46.008n, 


mjj/46.008 = 2n^ + n 2 


(C. 3) 
(C.4) 


Equilibrium relationship 


K 


K = 
P 


C’»o 2 , 2 ftojW 




(n 2 /n T ) P R 
(n 1 /n T ) 


\\ Ph 

_ ^ n 2 ^ P H 

n l n T 


(So ,) 2 P H 

V* 


' (C.5) 


(C. 6) 


The relationship defining pressure in terms of mass is derived as 
follows. Solving equation C.2 for n T 


n„ 


P V 
H H 

RT 


(C. 7) 


and combining with equation C.6 


_ n 2 RT 
K — T7 

P n ] _ V R 


(C.8) 


Solving equation C.8 for n^ 


n, = 


rt ; 

VK n 2 
H p 


Substituting C.9 into C.4 


(C. 9) 


"H 2RT 2 

n„ + n. 


46.008 V K “2 2 

H p 


(C.10) 


and solving C.10 for 


(* 


- 1 + 11+8 


n 2 “ 


Substituting C.ll and C.9 into C.2 

4RT 


RT 

V..K 46 
H P 


"h \ 1/2 


■h \ 
.008/ 


RT 

V H K n 

H p 


K 
_[ 
H 8 


P„ = -E- 


V„K 46.008 

H p 


0 


- 1+1 + 8 


RT 

V„K 46 
H p 


5.008 / 


(C.ll) 


(C. 12) 


A similar equation for the low pressure side may be written 


P L = 


K 

P 

4RT “l , / 

1 8RT 

\ \ i/2 i 

8 

V t K 46.008 \ 

L p ' 

1 + V t K 
^ L p 

46.008/ j 


(C. 13) 


Returning to equation C.l the following definitions can be made 


“ll " m H “ “ 


(C. 14) 


“L 




+ m 


(C . 15) 


where m^j is the total mass on the high pressure side at time zero 
and m^ is total mass on the low pressure side at time zero. 

m^ and m£ can be defined in terms of initial pressure using 
relationships defined in a manner similar to Equations C.12 and C.13. 
The differential Equation C.l can be written now as follows: 


Letting 


1H 

C 1L ' K P /8 ' 

(C. 16) 


4 RT 

(C. 17) 

2H = 

V„K 46.008 
H p 


4 RT 

(C. 18) 

2L = 

V t K 46.008 
L p 


C-3 


f = I P | C 1 h[ C 2H<^ - ”> - 1 Wl + 2C 2h (^ - 

- C 1L [ C 2L ( ”L + m) ' 1 + ' /l + 2C 2L ( "L + m) ] j (C ’ 19) 


This equation was solved numerically with various values of p using 
Milne's method to define m = m(t). A computer program incorporating 
the solution of the equation and the generalized least squares deter- 
mination of p was written for the analysis of the experimental data. 


